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ABSTRACT 

In this paper we consider the re-processing of high frequency photons emitted by He ii and 
He I during the epoch of cosmological recombination by He i and H i. We demonstrate that, 
in comparison to computations which neglect all feedback processes, the number of cosmo- 
logical recombination photons that are related to the presence of helium in the early Universe 
could be increased by ~ 40% - 70%. Our computations imply that per helium nucleus ~ 3 - 6 
additional photons could be produced. Therefore, a total of ~ 12-14 helium-related photons 
are emitted during cosmological recombination. This is an important addition to cosmological 
recombination spectrum which in the future may render it slightly easier to determine the pri- 
mordial abundance of helium using differential measurements of the CMB energy spectrum. 
Also, since these photons are the only witnesses of the feedback process at high redshift, 
observing them in principle offers a way to check our understanding of the recombination 
physics. Here most interestingly, the feedback of He ii photons on He i leads to the appearance 
of several additional, rather narrow spectral features in the He i recombination spectrum at 
low frequencies. Consequently, the signatures of helium-related features in the CMB spectral 
distortion due to cosmological recombination at some given frequency can exceed the aver- 
age level of ~ 17% several times. We find that in particular the bands around v ~ 10 GHz, 
~ 35 GHz, ~ 80 GHz, and ~ 200 GHz seem to be affected strongly. In addition, we computed 
the changes in the cosmological ionization history, finding that only the feedback of primary 
He I photons on the dynamics of He ii — > He i recombination has an effect, producing a change 
of ANe/Ns ~ +0.17% at z ~ 2300. This result seems to be ~ 2 - 3 times smaller than the 
one obtained in earlier computations for this process, however, the diff'erence will not be very 
important for the analysis of future CMB data. 

Key words: Cosmic Microwave Background: cosmological recombination, temperature 
anisotropics, radiative transfer 
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' 1 INTRODUCTION still be visible at mm, cm and dm wavelength. Observing these sig- 

^ , . natures from cosmological recombination may offer an indepen- 

It IS well known that cosmological recombination of hydrogen and , ^ . j , ■ r i i ■ i 

" r Qsnt way to determine some of the key cosmological parameters, 

helium in the Universe leads to the emission of several photons , , i- , , r- r , ^, r r ■ r 

such as the primordial helium abundance, the number density of 

per baryon, modifying the cosmic microwave background (CMB) , , , ^. tt, r , ■ ■ \r 

r — f= — T II — TTX I 7"; — II " I I ; — H baryons and the LMB monopole temperatu re at recombination (tor 

energy spectrum I Zeldovich et al. 1 96ol : Peebles 1 96a : DubrovicrJ TTTTTnrJTJT^ 

l /'l j f-'r— ^ . , n „ , 1 1 . J„ J, „ , . — , , overview e.g. see Sunyaev & Chluba^ 20091) . 

119751 : iDubrovich & Stolvarovl Il997h . Recently, detailed compu- ^ ' ^ ^ ^ 



tations of t he cosmological recombination spectrum were car- However, to approach this observationally challenging task it 

ried out (e. g.'Rubino-Ma rtm et al.l2006l : lchluba & Sunyaevl2006bl : is important to theoretically understand all the possible contribu- 
iRubifio-Ma rtin et al. 200^. showing that the recombinations of hy- tions to the cosmological recombination spectrum and how they 
drogen and helium lead to relatively narrow spectral features in may be affected by physical processes occurring at high redshift. 
the CMB energy spectrum. These features were created at redshift In this paper we investigate the re-processing of energetic pho- 
z ~ 1400, ~ 2500 and ~ 6000, corresponding to the times of Hi, tons initially released by Hei and Hee during cosmological re- 
Hei and Heii recombination, and, due to redshifting, today should combination by neutral hydrogen and helium. For this problem in 

particular quanta emitted in the Heii Lyman a (transition energy 
E ~ 40.8eV), Hei 2'Pi - I'So {E ~ 21.2eV), Hei 2'?! - I'So 
* E-mail: jchluba@cita.utoronto.ca (E ~ 21.0 eV) dipole transitions are important, but e.g. also the 
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higher Hei w'Pi - I'So series and Heii two-photon continuum do 
play some rolJl In contrast to the CMB spectral distortions created 
by helium at low frequencies, these photons lead to a large devia- 
tion of the CMB spectrum from the one of a pure blackbody, so 
that, after some (significant) redshifting, they are able to re-excite 
energetically lower-lying atomic transitions in H i and He i, signifi- 
cantly affecting the net rates of resonant and continuum transitions 
starting from the ground-state. This can lead to both changes in the 
cosmological recombination spectrum and the cosmological ion- 
ization history, and, as we demonstrate here, in particular the to- 
tal contribution of photons related to the presence of helium in the 
early Universe is increased by ~ 40% - 70% in comparison to com- 
putations that do not include the feedback processes under discus- 
sion here. Such a large addition to the cosmological recombination 
spectrum is very important, since in the future it may render a de- 
termination of the primordial helium abundance using differential 
measurements of the CMB energy spectrum slightly easier. 

To understand the physics behind this problem, we distinguish 
between two main types of feedback: (i) the self-feedback or intra- 
species feedback; and (ii) the inter-species feedback. The first type 
of feedback is related to photons that are emitted by some atomic 
species (e.g. Hei) and then affect the lower-lying transitions of the 
same atomic species. Since the difference between the time of emis- 
sion and feedback is connected with the redshift interval it takes 
to cross the energy distance between the resonances that are af- 
fected, for H I and He ii one therefore typically expects a delay of 
Az/z ~ 1% - 20%. In contrast to this, for the inter-species feedback 
energetic photons are emitted by some atomic species (e.g. He ii) 
and then feedback on lower-lying transitions of some other atomic 
species (e.g. Hei). Here the typical delay between emission and 
feedback is significantly larger (e.g. reaching Az/z ~ 40% in the 
case of Heii Lyman o- to Hei I'So continuum feedback as shown 
in Sect. |8.6l l, since the possible energy differences are much bigger. 

For h ydrogen the self-feedback problem was already st ud- 
ied earlier jChluba & SunvaevI |2007| : ISwitzer & Hiratall2008ah in 
connection with the Hi Lyman series. There, for example, pho- 
tons escaping from the Hi Lyman /? resonance, after redshifting 
by Az/z ~ 16%, will feedback on the Hi Lyman a line, leading to 
a small inhibition of hydrogen recombination. This occurs because 
the H I Lyman yS feedback adds photons to the phase space density 
around the Hi Lyman a resonance and hence increases the popu- 
lation of the 2p-state. These additional photons have to reach the 
very distant red wing of the H i Lyman a line (e.g. via redshifting 
or by some chain of transitions towards higher levels or the contin- 
uunj^that eventually leads to an emission in the 2s- Is two-photon 
channel) before the initial re-excitation of the electron by the feed- 
back is reversed. Similarly, photons emitted in the He i 2' Pi - 1 'So 
resonance will feedback on the Hei 2^Pi - I'So line, but here 
the difference between the emission and feedback redshift is only 
Az/z ~ 1%, due to the smaller energy distance between these reso- 
nanc es. The latter problem ha s as well been addressed in the litera- 
ture jSwitzer & Hiratj|2008 j) , also including the fact that on their 
way from one resonance to the other some of the He i photons are 
absorbed by the small amount of neutral hydrogen that is already 
present at redshift z ~ 2500 before they can actually feed back. 



' We give a detailed inventory of possible primary He i and He n feedback 
photons in Sect. l6.2l and l8.1l 

^ Excluding the Lyman a resonance itself, the main way out of the 2p- 
level is via a transition to the 3d state (occurring in ~ 90% of the cases), 
while a transition to the continuum happens with < 0.5% probability 
iChluba & Sunvaevll2009bl) 



In this context, it was also shown that the intra-species 
feedback of Hi and Hei photons leads to some small correc- 
tion to the cosmological ionization history, r eaching AN^/N^ ~ 
0.2% - 0.3% at z ~ 1100 for hydrogen jChluba & SunvaevI 
l2007h and ANJNc ~ -4% - 0.5% at z ~ 2000 for helium 
jSwitzer & Hiratal 120083) . This level of precision in our under- 
standing of the dynamics of cosmological recombination will 
be important for the analysis of upcoming CMB data from the 
Planck Surveyor, which was successfully launched in May this 
year. In particular, our ability to precisely measure the spectral 
index of CMB fluctuations may be compromised by the neglect 
of physical processes that can affect the ionization history at the 
level of ~ .1% close to the maximum of the Thomson visibil- 
ity function jSunvaev & Zeldoviclj|l970l) at z ~ 1100. Over the 
past few years many such processes have already been identi- 
fied ( for overview e.g. see iFendt et al.ll2009l: ISunvaev & Chlubal 



20091) by several independent groups (e.g. iDubrovich & Grachev 
2005; 'Chluba & Sunyaev 2006b; 'Kholupenko & Ivanchik* '200( 
Rubino-Martin et al. 2006; Wong & Scott 2007; Switzer & Hirat; 



2008b 



Karshenboim & Ivanovl 1200 8l ; iLabzowskv et ail 12009 



Jentschurdl2009h . also emphasizing that in principle all these pro- 



cesses do dir ectly change the cosmolo gical recombination spec- 
trum (e.g. see lChluba & Sunvaevir2009iJ) . 

Although the intra-species feedback affects the dynamics of 
recombination, it does not lead to any significant change in the 
cosmological recombination spectrum, and in particular the total 
number of (low-frequency) photons released during recombina- 
tion. This is mainly because the total number of available feedback 
photons is smalJJ but as we explain here (see Sect. 18. 3t . also the 
branching ratios to other levels play some important role. 

In terms of photon production the inter-species feedback is 
much more interesting, and was not taken into account at full depth 
so far. The simplest example is connected with the feedback of pri- 
mary He I 2^Pi - 1 ' So photons on H i. As we show here, practically 
all these photons never reach frequencies below the H i Lyman con- 
tinuum. Furthermore, we find that they feedback during the pre- 
recombinational epoch of H i. At that time the degree of ionization 
for hydrogen is still very close to the equilibrium Saba ionization 
(even when including the feedback). Obviously, the presence of ad- 
ditional non-equilibrium ionizing photons tends to increase the de- 
gree of ionization, but the recombination rate is high, so that this 
feedback eventually only leads to additional features in the cosmo- 
logical recombination radiation, but no significant correction to the 
ionization history. 

Here the most important aspect is that the electron which 
is liberated by the Hei feedback on the Hi Lyman continuum 
afterward has the possibility to recombine to some highly ex- 
cited state and then emit several photons on its way towards the 
ground-state. This is in stark contrast to the intra-species feed- 
back, where it is unlikely to reach very highly excited levels. The 
inter-species feedback therefore leads to loops of atomic transi- 
tions in the non-equilibrium ambient radiation field of the CMB 
jLvubarskv & SunvaevI [l983l) , which tend to erase the high fre- 
quency spectral distortion introduced by Hei. In these loops one 
energetic photon is absorbed, while several low-frequency photons 
can be emitted. A similar process was studied in connection with 
the release of energy by some decayi ng or annihilating particle s 
prior to cosmological recombination jChluba & SunvaevI l2009cl) . 



' e.g. th e number of H i Lyman /? relative to H i Lyman a photons is about 
~ 0.5% jChtuba & Sunvae\j|2007l) 
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As we explain here, per primary He i feedback photon about 2.6 y 
are produced by H i in addition to one H i Lyman a photon replac- 
ing the He i feedback photon (see Sect. l7.1.Tl for more explanation). 
We also study the feedback of Heii photons on both Hei and Hi, 
in detail explaining all the important physical aspects of this prob- 
lem (Sect.[8]l, and computing the effect of Heii Lyman a feedback 
(Sect. 18. 818. 9b on the cosmological recombination spectrum. Our 
main results for this problem are shown in Fig.l23landl24l 

The paper is structured as follows: in Sect. |2] we give the for- 
mulation and solution of several photons transfer equations which 
are important in the context of the He: intra-species feedback prob- 
lem. This Section is rather technical, and is only addressed to the 
interested reader. In Sect.[3]and|4]we apply the results of Sect.[2]to 
the problem of the He i 2' Pj - 1 ' So and He i 2^Pi - 1 ' So resonances 
including the H i Lyman continuum, solving them numerically. This 
Section provides some intuition for the important aspects of the 
problem, which then lead to the analytic approximations for the net 
rates and escape probabilities derived in Sect.|5] In Sect.|6]we com- 
pute the corrections to the ionization history and in Sect. I7I8I we 
discuss the changes in the cosmological recombination spectrum. 



2 DIFFERENT TRANSFER EQUATIONS IN THE NO 
SCATTERING APPROXIMATION AND THEIR 
SOLUTIONS 

In this Section we give the transfer equations and their solutions 
including different combinations of emission and absorption pro- 
cesses. In particular we are interested in the combined problem of 
the He I 2'P, - 1 ' So and He i 2^Pi - I ' So resonances and the H i Ly- 
man continuum. However, the equations and solutions given here 
can be easily applied to the other resonances of neutral helium and 
may also be useful for further studies related to the CMB spectral 
distortions generated by helium. 

We assume that the modifications to the solutions caused 
by partial frequency redistribution can be neglected, and work 
in the no line-scattering approximation. Thi s procedure has 
been u sed in several re c ent st u dies (e .g. ISwitzer & Hiratal 
l2008al : IChluba & SunvaevI l2008d. l2009dh and seems to be 
well justified given that the additional correct ions due to par- 
tial fr equency redistrib ution a re sm al l (e.g . 'Switzer & Hiratal 



2008al: [R ubiiio-M artm et alj l2008l : Iffirata & Forbe j |2009| ; 
Chluba & SunvaevI (2009a). However, as we comment here 



(see Sect. |6.4| l, when the considered resonances are very close to 
each other, additional coiTections are expected. 

We will then use the results of this section to compute the 
CMB spectral distortion due to the He i 2'Pi - 1 ' So and He i l^V^ - 
I ' So resonances at different stages of He ii He i recombination 
(Sect.|3}. This also allows us to check for additional corrections (i.e. 
due to time-dependence, the thermodynamic correction factor and 
line cross-talk including feedback) to the effective escape probabil- 
ity of these resonances with respect to the normal quasi-stationary 
approximation (Sect.|4ll. 



2.1 Isolated helium resonances 

Following the procedure described in IChluba & SunvaevI l l2009bh . 
in the no scattering approximation the transfer equation for photons 
in an isolated helium resonance connected with the transition from 



an excited level / to the He i ground state can be cast into the form 

1 dM; 
c df 



(la) 



Here is the one photon death probability of the resonance; 
c\ = -r- t¥- is the averaged one photon cross-section of the res- 

onance, where B^, is Einstein-B-coefficient of the line and Avp 
its Doppler width; N^^' is the population of the helium ground 
state; 4''.joJy,z) denotes the effective absorption profile; and N™'' 
describes the production of photons in the line. 

If we neglect two-photon corrections to the shapes of the pro- 
files in the different absorption channels contributing to the effec- 
tive absorption profile of the resonance, and if we assume that all 
are given by the normal Voigt profile 0y(v, z) for the considered 
resonance, then it is clear that 



fi - J_„h[v-v,ykTy 



(lb) 

(Ic) 



where /,', denotes the frequency-dependent thermodynamic correc- 
tion factor, which was discussed ear lier in connection with the h y- 
drogen Lyman a transfer problem jChluba & SunvaevI l2009dl fa). 
There it was shown that this factor follows from the detailed bal- 
ance principle and is required in order to fully conserve a blackbody 
spectrum in the distant wings of the line profile for the case of ther- 
modynamic equilibrium. Furthermore, one can directly write 

^iii^xi.^, (Id) 



where A'^^^^ is only redshift dependent. Here v, is the transition fre- 
quency of the resonance, gi, and g, are the statistical weights of the 
helium ground state and level i, respectively. In addition, and Rj 
denote the one photon rates at which electrons enter and exits the 
level ( via all possible channels excluding the resonance itself. 

If we now introduce the effective absorption optical depth in 
the resonance as 



Pd 



co-[N'^^' 



Jx[l+z],z)dz 



(2) 



where we defined the dimensionless frequency x 
we can write Eq. 0} as 



dA', 



v/(l + z), then 



(3) 



where f'^^_, = 9,7'^^^^. The solution of thi s equation in th e expanding 
Universe was already discussed earlier dChluba & Suiwaev 2009^). 
With the notation used here it can be written as 



AA';(z) = IKJz) - X Ff'(z). 



(4a) 



where the function Fy represents the frequency dependent part of 
the solution for the spectral distortion, which is given by 



Ff'iz) : 



&Az,z') = - 



J 0;,(z,z')'5-e-<bs<"~'~)dz' (4b) 
1 _ &[{z,z') 



(4c) 



where AA^„ = A',, - Nl\ V = x[\ -(- z'] and at z > Zs the CMB 
spectrum is assumed to be given by a pure blackbody spectrum Ny\ 
Furthermore, N'^Jz) = A';„(z)/[1 + zf, x, = v,/[I -H z], x; = v,/[I -H 
z'l and Nf = Ajr^^, with To = 2.725 K jFixsen & MatheJ 
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l2002h . Note that Nf does not explicitly depend on redshift. Also 
we have used that fy'{z')Nf = N'f, ■ 

For numerical computations it is important to analytically sep- 
arate the main contribution to the solution Eq. l l4bt . It can be ob- 
tained using the quasi-stationary assumption for the photon emis- 
sion rate (i.e. setting &[(z,z') = 1) and neglecting the thermody- 
namic correction factor in the definition of @'^(z,z'), which yields 



F^'-(z) = F-;\z) + ^F^/(z) 
1 _ e-4.<>--~' 



F^'iz) 



^Ff'iz) ■ 



J [0Uz,z)-l] d-.e 



dz', 



(5a) 
(5b) 

(5c) 



where @'.JyZ,z!) is defined by Eq. ( |4cl (. This solution describes the 
time-dependent evolution of the spectral distortion due to the emis- 
sion and absorption of photons in the resonance. The properties of 
this solution are very similar to the case of hydrogen, when neglect- 
ing two-photon corrections to the shape of the absorption profile 
dChluba & Sunvaevll2009bl) . We will discuss this solution in more 
detail in Sect. [3] 



2.2 Hydrogen Lyman continuum channel alone 

If we only consider the evolution of the photon distribution in the 
Is-continuum channel of the hydrogen atom then we can start from 
the equation describing the emission and absorption of photons due 
to direct recombinations to or ionizations from the ground state 
(e.g. see [Ch luba & Sunyaev 2006ai') 



1 dN^ 
c dt 



2v 



= N,N^MT,)a-i^--e'^ - <'o-i,,A',. (6) 



Here crisc is the Hi Is-photoionization cross section; N'^' is the 
ground state population of the hydrogen atom; the electron tem- 
perature, which is always very close to the radiation temperature 
Ty = Toil + z) with Tq = 2.725 K; A'e and A'p are the free electron 
and proton number densities; and from the Saha-relation one has 

(„Hi \LTE ^ 
1^) = {27^} ' ^'""^'^''^ Where hy^ ~~ 13.6 eV 
is the ionization energy of the hydrogen ground state. Note that in 
Eq. ^ we neglected the effect of stimulated recombination, since 
during recombination the photon occupation number around the 
Lyman continuum frequency is very small (n^ 1). 
Comparing with Eq. ( fTat and Eq {3} one can write 



df 



with the definitions 



N't! 



r 



■tri,,(x[l+z]X ' 
H{i+z) 



■Az 



2vl N,N^f,,{Ty) j_ 

\3/2 



(7a) 

(7b) 
(7c) 



/ 2 \ ' 

where fu(Ty) = L^J'^^^ 1 , and the thermodynamic factor = 
it.g/i[F-n.]/*ry jsjQjg (i^gj i^gj-g ijave assumed T^. = Ty. This as- 
sumption is very well justified at redshift z <; 1000 (e.g. see 
ISeager et alj|2000h . Now it is clear that the solution of Eq. Q is 
given by Eq. ([4j and similarly that Eq. l[5) is applicable. 



2.3 Cross-talk of the Hei 2'P| - I'So and Hei 2'Pi - I'So 



resonances 



If we now consider the problem for the simultaneous evolution of 
photons in the Hei 2' Pi - I'So and Hei 2''Pi - I'Sq resonances, 
then the transfer equation reads 



(8) 



where the superscript 'a' is related to the He: 2'Pi - I'So tran- 
sition and 'b' to the Hei 2^Pi - I'So transition. The solution of 
this equation can again be found using the same procedure as in 
IChluba & Sunvaevl l l2009j) . yielding 



AN^iz) = AA'^"(z) + AN^-\z) 
AN^-'iz) = IKJz) - K]] X F^^'iz). 



where the function F,, ' is given by 



Ff-^(Z) = £ @^(z,z:)^e-<Si- 



'dz'. 



(9a) 
(9b) 



(9c) 



with @l{z, z') defined by Eq. ^ and t^*^ = tI^^ + t|;^,. 

For numerical computations it is again better to analytically 
separate the main term in fJ' '(z). This results in 



F^'{z) = F[:'{z) + AF^--iz) 



AFf'iz) =J [0i,(z.z')e""* - l]a;-e-^*<"~~>dz'. 



(10a) 
(10b) 



where F[:'(z) is given by Eq. l lSbt and t^*^^' = t^^^ ■ 

Comparing Eq. l llObt with Eq. l l5c] l. one can see that the 
cross-talk between the lines simply leads to a frequency- and 
time-dependent modulation of the effective photon emission rate 
0L(z, z') related to the resonance ; by Physically this 

just reflects the fact that with time some photons from the resonance 
/ will disappear from the photon distribution due to the absorption 
in the other resonance. 



2.4 Hei 2'Pi - I'So and Hei 2'Pi - I'So resonances with 



hydrogen continuum opacity 



With the results of Sect. l2.3l It is now straightforward to write down 
the solution for transfer problem in one helium resonance including 
the hydrogen continuum. This yields 



AA'::'^(z) = AA'^=-"/''(z) + AiV^«(z) 
AA?-''"(z) = [N[^^(z) - <] X Ff{z). 

where the function F^y'' can be written as 



F["(z) = f;;'(z) + AF["(z) 



(11a) 
(lib) 

(12a) 



AF--'(z) =J~[0;,(z,z')e"^"^'' 



^'^-V.;'.;) _ I \q_, e-^s*"-'-) dz', (12b) 



T r + t"-, - t', . Here 

abs abs abs 



with 0i(z,z') defined by Eq. (gcli; rj^*"'' - 
we used the notation 'a/b' in the superscripts, which means either 
resonance 'a' or 'b'; and F[;'(z) as defined in Eq. l lSbt . 

Furthermore, when simultaneously including both resonances 
and the H i continuum one finds 

ANf (z) = AA^f -"(z) + AA^f (z) + AN'f-'' (z) (13a) 



AA'f ''(z) = [K^{z) - <'] X Ff ''(z). 



(13b) 
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where the function F^" ' is given by 

Ff-\z) = F•;\z) + ^Ff■^{z) (14a) 
KFf'\z) =^\®'Sz, zOe""*""*"" ~* - l]a,-e-^*("~ .=»dz', (14b) 

with ©- (z, z') defined by Eq. r-""- = r»,^ + r^,^ + r^,^ - r^,^; 
and Fy^iz) as defined in Eq. dSbt 

Note that it is very easy to extend the solution to n differ- 
ent resonances with Hi continuum absorption between the lines, 
or other processes of emission and absorption (e.g. Hi 2s-ls two- 
photon emission). However, since in most cases resonances in hy- 
drogen and helium are very distant to each other, for our purpose 
we neglect the case of cross talk between many resonances. 



3 NUMERICAL RESULTS FOR THE HIGH 

FREQUENCY CMB SPECTRAL DISTORTION AT 
DIFFERENT TIMES 

Using the solutions given in the previous sections it is possible to 
compute the high frequency spectral distortions due to the emis- 
sion of photons in the main helium resonances at different epochs. 
To demonstrate the main aspects of the problem here we shall only 
consider the evolution of the photon field due to the He i 2 ' P i - 1 ' So 
and He I 2^^Pi - 1 ' So transitions. These two resonances are separated 
by Av/v ~ 1%, so that photons from the Hei 2'Pi - I'So line will 
mainly feedback on the Hei 2'Pi - I'So intercombination line af- 
ter Az/z ~ 1% which corresponds to about 700 Doppler width at 
z ~ 2200 (see Table [T]!. In addition, it is expected that due to the 
presence of a small fraction of neutral hydrogen during helium re- 
combination resonance_ghotons_wnb Lyma n 
continuum dSwitzer & Hiratal l2008d : | Rubino-Martm et ai] |2008h . 
These pho tons will re-appear as pre-r ecombinational photons from 
hydrogen jRubiiio-Martm et"al]|2008l) . at frequencies very far red- 
ward of the main helium resonances. This leads to two effects: 

(i) At a given time the main resonances will be supported by 
fewer photons than in the case without H i continuum opacity. This 
leads to an increase in the effective esca pe probability which be- 
comes very large at redshift z £ 2400 iKholupenko et al.l [20071 ; 
ISwitzer & Hiratall2008 jjRubiiio-Martm et alj2008l) ^ 

(ii) The feedback of photons from the Hei 2' Pi - I'So resonance 
on the Hei 2^Pi - I'So intercombination line will be reduced, 
since part of the photons released in the Hei 2'P| - I'So line 
will disappear from the photon distribution before they can actu- 
ally reach the Hei 2^P| - 1 ' So resonance (see also explanation in 
ISwitzer & Hiratall2008al) . As we demonstrate here due to this pro- 
cess the feedback of /(Hei 2'Pi - I'So) ^ r(HEi 2'Pi - I'So) 
practically stops at redshift z g 2200 - 2400 (see Sect. 16. 3t . 

In order to compute the solution for the spectral distortion 
caused by the considered resonances one has to give the solutions 
for the population of the He i levels as a function of time. For this 
one has to assume that the approximations used in the multi-level 
helium recombination codqj already captures the main processes, 
and that the corrections due to the additional effects (e.g. time- 
dependence, thermodynamic correction factor, line feedback and 
line cross-talk) are small. 



It was already shown i Kholupenko et alj 1 20071 : 
ISwitzer & Hiratj l2008al : iRubifio-Martm et alj |2008|) that the 
main correction during helium recombination in comparison 
with the standard computation dSeager et al.l |2000|) is due to 
the speed-up of the Hei 2'Pi - I'So and Hei 23p, - I'So 
channels caused by the Hi continuum opacity. To include this 
process in our computa tions of the Hei po p ulatio ns we will 
follow the approach of iRubino-Martm et ak I ( l2008l) using the 
ID-integral approximation (see Eq. (B.3) in their paper) to take 
the increase in the escape probability of the Hei 2'Pi - I'So 
and Hei 2'Pi - I'So transitions into account. In addition, for 
the Hei 2' Pi - I'So resonance, corre ctions due to par t ial fre- 
quency redistribution are important dSwitzer & Hiratj l2008d: 
iRubifio-Martm et alj2008l) which we accoun t for with the ' fudge '- 
function used in iRubifio-Martm et"aL I l l2008h . Henceforth we will 
refer to this mo del as our reference model. Consistent with our 
previous works jRubino-Martm et alj l200d: IChluba et al.l |2007| : 



iRubifio-Martfn et alj|2008l : IChluba & Sunvaevll2009bl lal) we used 
the cosmological parameters Tq = 2.725 K, Fp = 0.24, h = 0.71, 
fib = 0.0444, fi,„ = 0.2678, Ha = 0.7322, and Q.^ = 0. 

To understand the importance of the different corrections 
for the shape of the CMB spectral distortion introduced by the 
He I 2' Pi - 1 ' So and He i 2^Pi - 1 ' So transitions we now discuss the 
solution for two different representative stages (Sect. [3TT] and |3.2| > 
during helium recombination. We then also compute the present- 
day (z = 0) CMB distortion at high frequencies including differ- 
ent processes. In particular we show that the absorption in the Hi 
continuum completely erases the high frequency spectral distortion 
from H e n — » He i recombination , a fact that was already sus pected 
earlier l lChluba & Sunvaevll 20071 : iRubiiio-Martm et alj200^ . 



3.1 High frequency spectral distortion at z = 2400 

In Fig.[T]we show the results for F,, at z = 2400 for different combi- 
nations of the line emission and absorption processes. We normal- 
ized all the curves to the central distortion in the Hei 2' Pi - I'So 
line, AA';|_(z) = AA'^^,(z) = N%^{z) - N^^. In all cases we neglected 
the possible emission in the H i continuum channel, since close to 
the resonances it contributes very little to the total distortion. For 
comparison we also show the simple quasi-stationary solution 



FT-'-'\z) 



(15) 



^ We based this code on the works of iRubiiio-Martm et all j200l 
IChluba et all )2007l) . and lRubiiio-Martm et alj )2008h As s imilar code w 
already used for the first training of Rico lFendtetalJl2009h . 



which neglects any time-dependence, cross talk, the H i continuum 
opacity, or feedback of the lines. Here r|j = Tj, where Tj is the 
standard Sobolev optical depth for the resonance ;, and iS.N\^{z) = 
ANy^(z) is the central distortion in the Hei 2^^ Pi - I'Sq resonance. 
Also we defined^',, = ip'(v')dv' with (/iy(v, z) = 0v(^'Z)/Avp. 

Note that the effective absorption optical depth, t[|, takes into 
account that only a fraction p'^ of interactions with the resonance re- 
ally leads to an absorption or death of the photon. As an example, 
for the H i Lyman a resonance this absorption is related to a tran- 
sition of the electron towards higher levels or the continuum in a 
two-photon process (e.g. see Ichluba & Sunvaev 20091j), however, 
for the higher Hi Lyman-series also spontaneous decays towards 
lower levels matter All these contributions can be taken into ac- 
count using the appropriate branching ratios for each excited level 
in hydrogen and helium, assuming that the ambient radiation field 
is given by the CMB blackbody spectrum. 

When considering the two resonances separately and neglect- 
ing the H I continuum opacity (Fig.[T] upper left panel), F,, is simply 
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Figure 1. Spectral distortion due to tlie Hei 2'Pi - I'So and Hei - I'So resonances at z = 2400 - upper left: isolated lines - upper right: cross-talk 
without H I continuum opacity - lower left: cross-talk with H i continuum opacity - lower right: comparison of the total distortion. The frequency is given 
as the distance to the central frequency of the Hei 2'Pi - I'So resonance in Doppler units, xq = ^rr-- In all panels we nonnalized the total distortion 

ANy{z) by AN'^^iz) = A'em(z) "^nl for the Hei 2'Pi - I'So resonance, i.e. Fy(z) = ANy{z)/ AN'^^{z). Furthermore, Fy(z) denotes the total distortion, F'fXz) the 
distortion due to the Hei 2'Pi — 1 'So line, and Fy(z) the distortion due to the Hei 2^Pi - 1 ' So resonance. For comparison we show the simple quasi-stationary 
solution FJJ''''*(z) according to Eq. il5\ . 



given by the sum of the distortions from each resonance. At low fre- 
quencies one can clearly see the modulation of the distortion due to 
the time-dependenc e of the emission rate, like i n the case of the Hi 
Lyman a distortion jChluba & Sunvaevll2009d) . In addition for the 
distortion caused by the Hei 2'Pi - I'So resonance one can also 
find the scaling Ff, ~ i/fy, aga in in full analog to the H i Lyman a 
line dChluba & Sunvaevli2009bl) . 

If we now include the cross-talk between the lines (Fig. [T] 
upper right panel), then the main difference is only due to the 
fact that the photons from the Hei 2' Pi - I'So line have to pass 
through the He: 2^Pi - I'So resonance, when they feedback on 
the Hei 2'Pi - I'So transition and get partially re-absorbed. This 
leads to a drop of by £""^11 in the vicinity of the He i 2' Pi - 1 ' So 
resonance (at around .^d 700). However, the additional absorp- 



tion of photons from resonance ; in the (distant) damping wing of 
the line j is negligible. Since tJ] is a function of redshift, the am- 
plitude of this drop depends on the redshift of line crossing. The 
farther one goes red- ward of the Hei 2^ Pi - I'So line the smaller 
the absorption caused by the line crossing becomes. At xo < -10"* 
the distortion Fj) again becomes comparable to in the case with- 
out cross-talk (compare curves in the upper panels of Fig.[T]l. This 
is because those photons have passed through the Hei 2^^ Pi - I'So 
resonance at much earlier times, when the effective He i 2^Pi - 1 ' So 
absorption optical depth was smaller (cf. Fig. |2j. At z = 2400 for 
our reference model one has rjj ~ 0.8, so that the distortion due 
to the Hei 2'Pi - I'So line should drops by a factor of ~ 2.2 at 

Xd 700. This is in good agreement with our computations (cf. 

Fig-El upper right panel). 

Finally, when we also include the effect of the Hi contin- 
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Figure 2. Effective absoiption optical depth, = PjTj, in the first two 
Hei ji'Pi - I'So intercombination and He: n'D2 — I'So quadrupole tran- 
sitions. In all cases the death probabilities p'^ are very close to unity, so that 
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Figure 3. Absorption optical depth, r^, in the Hi Lyman-continuum. All 
curves can be computed using the definition Eq. )7bt with the Hi Is- 
photoionization cross section, where the appropriated Az/z. is defined by 
the distance between the considered resonances. The dotted curve shows 
the optically depth when crossing the Doppler core of the Hei 2'Pi - I'So 
resonance. The other curves show the total optical depth between the quoted 
resonances. Thes e have been computed usin g the solution for N^^ from our 
reference model iRubiiio-Martin et alj20O8h using Eq. )7bt . 



uum (Fig. [T] lower panels), the shape of the distortion changes 
drastically. In particular at large distances red-ward of the reso- 
nance the remaining distortion due to Hi absorption is strongly 
reduced. From this Figure it is clear that basically all the pho- 
tons from the Hei 2'Pi - I'So and Hei 2^Pi - I'So lines are re- 
absorbed in the H i continuum before they can actually reach fre- 
quencies below the Hi Lyman continuum threshold frequency v^, 

which at z ~ 2400 is at xd 2.1 x 10^. These absorbed photons 

should lead to additional ionizations of neutral hydrogen well be- 
fore the actual epoch of hydrogen recombination (z < 1600). Since 
at those times hydrogen is still in very close equilibrium with the 
continuum, this process will not cause any important changes in 
the ionization history, but should lead to some pre-recombinational 
emission in the lines of hydrogen. However, until now this effect 
has not been fully taken into account, since only those absorp- 
tions causing the changes in the effective escape probabilities of the 
Hei 2'Pi - I'Sp and Hei 2'^Pi - I'So resonances where accounted 
for ( iRubirio-Martfn et"al]|2008l) . We will discuss this problem in 
more detail below (Sect.|7]l. 

Again looking at the lower left panel in Fig. [T] one can also 
observe some modifications of the spectral distortion between the 
Hei 2' Pi - 1 'So and Hei 2^Pi - I'So resonances. These changes are 
cause by the absorption of photons from the Hei 2' Pi - I'So line 
while they are on their way to the He i 2^Pi - 1 ' So intercombination 
line. This will lead to a reduction of the feedback correction to the 
Hei 2^Pi - I'So transition, as we will explain in detail below (see 
Sect. l6.3t . One can estimate this reduction by simply computing the 
H I absorption optical depth between the two lines (see solid line in 
Fig.O. At z ~ 2400 we find t"^'' ~ 0.7, so that, in good agreement 
with the result shown in Fig.[T] until the Hei 2'P| - I'So photons 
have reached the Hei 2"^Pi - I'So resonance due to Hi continuum 
absorption one expects an reduction of the spectral distortion by a 
factor of ~ 2. 



3.2 High frequency spectral distortion at z = 2000 

Moving to z ~ 2000 one can see that the general behavior of the 
solution is very similar to the case we considered in the previous 
section. However, for example, now the drop in the spectral distor- 
tion from Hei 2'Pi - I'So line the due to the feedback absorption 
in the Hei 2^Pi - I'So resonance has increased to about a factor of 
~ 8.5 (cf. Fig.|4] upper right panel). Again looking at Fig.[2]we find 
t[J ~ 2.1 at z = 2000, which confirms this result. 

Furthermore we can see that the influence of the hydrogen 
continuum opacity on the shape of the spectral distortion has be- 
come very drastic (Fig. |4] lower panels). The photon distribu- 
tion is practically narrowed down to those photons appearing in 
the close vicinity of the two resonances. Practically none of the 
photons emitted in the Hei 2' Pi - I'So transition really reach 
the Hei 2'Pi - I'So intercombination line. This implies that the 
feedback correction to the Hei 2"^Pi - I'Sq transition is already 
expected to be negligible. In Sect. 16.31 we will show that the 
Hei 2'Pi - I'So ^ Hei l^'Pi - I'So feedback practically stops 
at z < 2000 - 2200. 

Given that photons supporting the flow of electrons to the 
Hei 2' Pi and Hei 2'Pi via the considered resonances are only 
present in a narrow range of frequencies around v ~ v.^ and v ~ vj, 
it is also clear that time-dependent corrections and correction due 
to the thermodynamic factor /J cannot be very important. This is 
in stark contrast to the Lyman a escape problem during hydro- 
gen recombination, where a very large part of the total correc- 
tion to the effective escape p robability is caused by these processes 
dChluba & Sunvaevll2009bl) . This is because during hydrogen re- 
combination the Hi Lyman a spectral distortion is always very 
broad, so that even photons at |xdI ~ 10' - lO** matter at the re- 
quired level of precision. However, during He 11 — > Hei recombi- 
nation these contributions turn out to be negligible (see Fig. [6] and 
Sect.H. 
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Figure 4. Same as 

In Fig.|4]we also show the total spectral distortion at frequen- 
cies below the threshold of the Hi Lyman continuum. It is clear 
that only a very small amount of photons really reach below this 
frequency. In a complete treatment one should allow these addi- 
tional photons to feedback on the transition in hydrogen during 
the pre-recombinational epoch. However, the total amount of pho- 
tons created from this feedback is rather small in comparison to the 
emission produced from the absorption of helium photons before 
they pass the H i Lyman continuum. 



3.3 CMB spectral distortion at redshift z = 

In this Section we show the spectral distortion of the CMB caused 
by the emission and absorption in the Hei 2'Pi - I'So and 
He I 2'Pi - I'So lines as it would be visible to today for different 
assumptions on the considered processes. However, it turns out that 
when including the absorption in the H i Lyman continuum practi- 
cally all the photons released by these transitions are erased before 
they can reach below the threshold frequency of the H i Lyman con- 



10 shell Hydrogen and 10 shell Helium atoms, ; = 2000 
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ig.EJbut at; = 2000. 

tinuum. As we will explain in more detail below (Sect.|7}, all these 
absorbed photons from helium lead to additional emission by hy- 
drogen during its pre-reco mbinational epoch, which w as neglected 
in previous computations jRubiiio-Martm et alj2008l) . 

In Figure |5] we give the present-day (z = 0) CMB spec- 
tral distortion Aly caused by the Hei 2' Pi - I'So (left panel) 
and Hei 2^Pi - I'So (right panel) transition. In the computations 
we used the results for the populations from our reference model 
jRubino-Martm et al. including 10 shells in hydrogen and 

10 shells in helium. The dash-dash-dotted lines give the distor- 
tion computed in the normal 5-function approximation for the line- 
profile, when neglecting the speed-up of the helium recombination 
dynamics by the presence of neutral hydrogen. When also includ- 
ing this pro cess we obtain the solid cu rves, which were already 
presented in lRubino-Martm et alj ( |2008|) . One can clearly see how 
due to the speed-up of He i recombination caused by H i absorption 
in the latter case the high frequency spectral distortion from the 
considered lin es narrows significantly in comparison to the former 
case (compare lRubifio-Martfn et alj|2008l) . 

If we now compute the distortion including the time- 
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Figure 5. CMB spectral distortion due to the He i 2^Pi - 1 ^Sq (left panel) and Hei 2"^Pi - 1 ^Sq (righ t panel) transitions at ^ = f or different cases (see Sect. |3.3| 
for details). In all computations we used the solution for the population from our reference model l lRubifio-Martm et al.||2QQ8l) . However, for comparison we 
also considered the case without the speed-up of helium recombination caused by the hydrogen continuum opacity (dash-dash-dotted line). 



dependence of the emission process and the thermodynamic cor- 
rection factor, we obtain the dashed curves. For the He i 2^ Pi - 1 ' So 
resonance this leads to an increase of the spectral distortion in the 
frequency range 1800 GHz < v < 3200 GHz by about 10%, while 
the He I 2-^Pi - I'So line is hardly affected. Like in the case of 
the H I Lyman a line the wing contributions to the escaping pho- 
ton distribution for the Hei 2^ Pi - l^So are important, since in 
the Hei 2^Pi - l^So Doppler core the effective absorption optical 
depth is ext remely large, so that hardl y any photon emitted there 
can survive dChluba & Sunvaevll2009iJ) . Therefore in particular the 
scahng of the thermodynamic correction factor at frequencies be- 
low the line center ( f^ < 1) leads to additional leakage of pho- 
tons, increasing the effective escape probability . This is completely 
analog to the case of the Hi Lyman a line jChluba & SunvaevI 
l2009dl lbl). In contrast to this, the wing contributions in the case of 
the Hei - I'So resonance are not very important, since the 
Hei 2^Pi - I'So Doppler core only becomes mildly optically thick 
to absorption during helium recombination (cf. Fig. O. Therefore 
neither time-dependence nor t 1 can affect the escaping number 
of photons and hence the amplitude of the CMB spectral distortion 
very much. 

If we now in addition allow for cross-talk between the lines, 
we can see that the feedback absorption during the passage of 
Hei 2'Pi - I'S,, photons through the Hei 2^Pi - 1' resonance, 
as expected, leads to an e^'^J suppression of the spectral distor- 
tion, which becomes more important toward higher frequencies, i.e. 
lower redshift of photon emission. In the case of the He i 2'Pi - 1 ' So 
line one can also see a very small modification in the amplitude of 
the total distortion. This is due to the small amount of re-absorption 
of photons in the distant red damping wing of the Hei 2'Pi - I'So 
resonance, however this modification is negligible. Note that here 
we have not yet included the change in the dynamics of helium re- 
combination caused by the feedback of He i 2' Pi - 1 ' So photons in 
the Hei 2^Pi - I'Sq line, but the final correction is small, so that 
the correction to correction can be neglected. We will discuss this 
case in more detail below (Sect. |631 ). 

If we finally also include the H i continuum opacity in the com- 
putation of the spectral distortion, we can see that basically all pho- 
tons disappear. Only a very small number of photons emitted a very 



early times during helium recombination can escape. This addi- 
tional huge reduction of the spectral distortion due to H i absorption 
was not yet taken into account, and it leads to a pre-recombinational 
feedback to hydrogen, which induces the emission of several addi- 
tional H I photons during helium recombination, as we will show in 
more detail below (Sect.|7}. 



4 CORRECTIONS TO THE NET RATES AND ESCAPE 
PROBABILITIES: NUMERICAL RESULTS 

We now want to include the corrections in the solutions of the ra- 
diative transfer problem into the computations of the helium and 
hydrogen recombination history. Knowing the solution, ANy(z), of 
the CMB spectral distortion as a function of time, and assuming 
that the differences in the net rates of the lines introduced by the 
different corrections are small, one can incorporate the effect of 
these processes into the multi-level code by modifying the escape 
probabilities for the He i resonances and the H i continuunjf[ Here 
we will give the results of our numerical computations using the 
solutions of the photon transfer equation given in Sect.[2l In Sect.|5] 
we will give simple analytic expressions which then allow to in- 
corporate these corrections into the multi-level recombination code 
with sufficient precision. 



4.1 Corrections to the escape probabilities of the main 
helium resonances 

If we start with the helium resonances and neglect two-photon cor- 
rections as above then it is clear that the net change in the popula- 
tion of level ; by transitions to the ground state can be cast into the 



^ This is only one possible approach, which is equivalent to directly includ- 
ing the relative changes in the net rates of the lines. 
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Figure 6. Effective escape probabilities in the Hei 2'Pi - I^Sq (left panel) and Hei 2^Fi - I^Sq (right panel) resonance. We computed the effective escape 
probabilities using formula ill) and the solution for the helium and hydrogen populations from our reference model. For comparison we also give the standard 
Sobolev escape probabilty and the result of different analytic approximations (s ee Sect.|5]for more details). In the left panel we also show the relative difference 
of our numerical result with the ID-approximation jRubiiio-Martfn et al.l2008h as inlay. 



form (cf. lChluba & Sunvaevll2009hl) 



IT 



ridNy 

J c dt 



dvdO. 



Eq. [13 

I 



-p; hv.Bi,, N»:' {k,, - <} X P-^^. (16) 



Here the effective escape probability is defined as 



= j| v'v(y,z) [l - f;Xz)G[(z)] dv 



(17a) 
(17b) 



where AA'^ is the overall spectral distortion of the CMB at red- 
shift z introduced by all the considered line emission and absorp- 
tion processes. Furthermore, we have AN'^^^z) = A'cmfe) ~ ^v\z), 
and ip'yiv,z) = 0y(v, z)/Avq. Note that P'^giz) can be directly com- 
pared to the standard Sobolev escape probability P'^ using the re- 
latioti P[„Az) = p:P[M)/\1 - PlPUz)] with = I- p\ 
(e.g. see lChluba & Sunvaevll2009cl) . Below we will give the results 
for P'^giz) when including different combinations of the considered 
modifications. Note that for most practical purposes one can use 
^eff.sfe) ~ P'i Kff(z)' since usually pi^ P^(z) « 1. 

4.1.1 Numerical results for the Hei 2'Pi - I'So resonance 

First we consider the escape of photons from the Hei 2' Pi - I'So 
resonance. From the results of Sect. [3] it is already clear that one 
does not expect very large corrections due to time-dependent as- 
pects of the problem or the thermodynamic corrections factor. This 
is because the H i continuum opacity (even at rather early stages of 
helium recombination) has such a large impact on the shape of the 
spectral distortion (cf. Fig.lTJ. In addition it is clear that the cross- 
talk with the He: 2^^Pi - I'So resonance cannot be very important, 

* This relation is obtained using the quasi-stationary solution for the popu- 
lation of the initial level. 



since in comparison with the line center the Hei 2' Pi - 1 'So absorp- 
tion profile drops by several orders of magnitude until 700, 

and since there is hardly any photons from the Hei 2^Pi - I'Sq 
line close to the He: 2'Pi - I'So resonance. The feedback of pho- 
tons emitted by higher transitions in helium on the blue side of the 
Hei 2' Pi - I'So resonance will be much more important. We will 
account for this effect in an approximate way in Sect. 16.41 

In Fig.|6]we present the numerical results of our computation. 
In the case of the Hei 2' P, - I'So the relative difference to the ID- 
integral approximation for the modi fied escape probability wh en 
including the H i continuum opacity l lRubino-Martm et"aLll2008l) is 
only of the order of a few percent, where part of the correction is 
due to the ID-integral approximation itself. For the dynamics of 
helium recombination the impact of such correction is negligible. 
Here we will not consider them any further, however, there is no 
principle difficulty to take them into account. 



4.1.2 Numerical results for the Hei 2'Pi - I'So resonance 

As mentioned above (Sect.[3]l for the Hei 2'Pi - I'So transition 
corrections due to time-dependence and the thermodynamic cor- 
rection factor are not important. However, the feedback of photons 
from the Hei 2' Pi - I'So resonance should have a strong impact, 
in particular at early stages of helium recombination, when the H i 
continuum opacity between the two resonances is not yet becoming 
too large, so that most of the Hei 2'Pi - I'So photons do reach the 
Hei 2'Pi - I'So line center. 

In Fig.[6]we present the numerical results of our computation 
for the escape probability in the Hei 2^P| - I'Sq. One can see 
that when including the effect of the hydrogen continuum opacity 
but neglecting the cross-talk between the lines, the escape prob- 
ability closely follows the c urve computed with the I D-integral 
approximation given earlier ( iRubino-Martfn et"al]|2008l) . There is 
some 1% - 4% difference between the analytic approximation (di- 
amonds) and the numerical result for this case (crosses) at redshift 
z ~ 1800 - 2200, which could be avoided when using the full 2D- 
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integral expression given by iRubirio-Martm et alj ( |2008|) (Eq. B.l 
in their paper), but again sucii precision is not required for the re- 
combination dynamics of helium. 

When also taking the cross-talk between the lines into ac- 
count, the effective escape probability is strongly reduced at z > 
2200. As we will see below, this is mainly due to the feedback 
of He: 2'Pi - I'So photons on the Hei 2'Pi - I'So transition. 
Also one can see, that at z < 2200 the effective escape probabil- 
ity again follows the case of no feedback. As we will explain in 
Sect. |6.3| this is due to the fact that because of the strong absorption 
of Hei 2'Pi - I'So photons in the Hi continuum the feedback is 
simply stopped. For comparison we also show the simple analytic 
approximation, Eq. ([33}, which includes the effect of feedback from 
the He I 2'Pi - 1 ' So line and the absorption of photons between the 
line (dash-dotted curve). Again the agreement is sufficient at the 
desired level of accuracy. For a derivation of this approximation we 
refer the reader to Sect.O 



4.2 Escape and feedback in the H i Lyman continuum 

If we now consider the escape problem in the H i continuum, then 
we can define the net rate connecting the Is-state with the contin- 
uum as (see Appendix iBlfor derivation) 



"dT 



1 dA',, 
c "dT 



dvda 



(18) 



where and iJj'^ are the photo-recombination and -ionization 
rates of the Is-state in ablackbody ambient radiation field. The con- 
tinuum escape probability is also given by the expressions Eq. l ll7b . 
with appropriate replacements (see derivation in Appendix |B). In 
particular one has to replace (^^(v, ;) with 



(v) 



(19) 



5 CORRECTIONS TO THE NET RATES AND ESCAPE 
PROBABILITIES: ANALYTIC CONSIDERATIONS 

In this Section we now discuss the analytic approximations that 
can be used to include the feedback process into the multi-level 
recombination code. As we already argued in Sect. |4l for the he- 
lium recombination problem time-dependent corrections are not so 
important. Here we therefore use the the quasi-stationary approxi- 
mation. If necessary it is straightforward to compute more refined 
cases using the analytic expression given above. However, for more 
detailed computations also other corrections, (e.g. related to partial 
frequency redistribution and electron scattering) will probably be- 
come more important (see discussion below). 

We will consider to types of feedback, first the feedback be- 
tween difi'erent resonance of neutral helium. This will lead to a de- 
lay of helium recombination, which due to the additional absorp- 
tion in the H i Lyman continuum is sup pressed very st rongly. This 
problem was already discussed earlier in lSwitzer & Hirata (,2008ai) , 
however our final correction to the ionization history seems to be 
smaller (see Fig.llOll. 

The second type feedback is the one on the H i Lyman con- 
tinuum which leads to additional pre-recombinational emission 
from hydrogen during the epoch of He ii — > He i recombination. 
Part of this correction was already taken into account earlier by 
iRubiiio-Martfn et al] ( |2008|) . but as we have seen in Sect. [33] (e.g. 
Fig. [5] solid lines) in their computations some photons still escaped 
until z = 0. All these photons will still be re-absorbed in the Hi 
continuurrQ and lead to additional feedback, but as we show below 
(Sect. |63t , no net change in the ionization history. 



5.1 Net rates for tlie multi-level code 

For our multi-level code the net rates in the resonances are im- 
portant. Here we are interested in those resonances leading to the 
ground state of hydrogen or helium. The change in the popula- 
tion A^, of level / due to the transition i Is is given by (cf. 
IChluba & Sunvae v 2009dt) 



dA^, 



gH 



(20a) 
(20b) 



Note that ip'J^v, z) is normalized like f'Jy'^' z) dv = 1 . 

With equation ( 117b it now in principle is possible to compute 
the escape probability for the H i Lyman continuum including pos- 
sible time-dependent corrections. However, it is already clear that 
the direct escape in the Hi Lyman continuum is very small, and 
that it in particular does not lea d to any large correction to the re- 
combination dynamics (e.g. see lChluba & Sunvaevll2007l) . For the 
direct escape of photons it is therefore s ufficient to use the ana- 
lytic expression dChluba & Sunvaevll2007h P^^. « I/[l -I- rf"] with 
.j-csc _ kj^ estimates in our computations. A similar 

expression can be applied to include the direct escape of photons 
in the Hei I'Si continuum, but again this has a very small impact 
on the dynamics of helium recombination. In addition, for detailed 
computations one should simultaneously includ e the effect of the 
hydrog en continuum opacity into this problem dSwitzer & Hiratal 
l2008al) . 

However, the feedback from photons that were emitted by the 
main resonances of helium still has to be taken into account and 
does lead to some pre-recombinational emission by hydrogen, as 
we will explain in more detail below (Sect. [6^ 



where A,is denotes the Einstein A-coefficient of the transition / — > 
Is, h' = Y^N' is the average occupation number of the photon field 
over the line absorption profile. Furthermore, we have introduced 
the line occupation number n', = -^N', = - 1] » ^^Ht^. 

The last approximation is possible when induced effects can be 
neglected, however, in particular during the pre-recombinational 
epochs of the considered atomic species one should use the full 
expression for n[^. Note that for Eq. i20\ it is still assumed that the 
effect of stimulated transitions can be captured using equilibrium 
values (i.e. [1 + n'^] a [1 + n^']). 

In the Sobol ev approximation for the is olated resonance one 
will have (e.g. see lChluba & Sunvaevll2009dl) 



(21) 



where AAf[ = A'^ - K' and P'^ 



is the Sobolev escape prob- 
ability. On the other hand, in the no line scattering approximation 



This was already pointed out bv lRubifio-Martm et alj j2008h . 
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one finds (e.g. see IChluba & Sunvaevll2009dh 



em 



(22) 



The expression l l25b then implies 

AA'4„(z) 



rill-Mi] 



where = 1 - pj, and P'^ 



with r'. 



Tg. Note that in 



and with Eq. i fTTt for /,:,' = 1 



the limit p'^ —> I one has iV^.,. = N'^^^. It is valid for the low proba- 
bility Hei n^Pi - I'So intercombination lines and Hei ?7'D2 - I'Sq 
series during helium recombination. Physically this limit is equiv- 
alent to the approximation of complete redistribution in the line for 
each scattering event. However, here the complete redistribution is 
achieved via transitions to higher levels rather than attributing it to 
a normal scattering (i.e. Is — ^ 2p — » Is) event. 

Also one should point out that for A t[j » 1 one again finds 



Pi 



AN'Jz) 



(27) 



(28) 



so that ANJ(z) = NJ(z) - K] = AN^Jz) x [1 - P^]. 

To include this into our multi-level code, we want to replace 
AA'em(<:) with the expression related to AA'j^ = N[- Ny'.. For this we 
assume that the solution for the population of the considered initial 
leve l / is given by the quasi-st ationary solution. This implies that 



N' 

eailier in the case of the no line scattering approximation. Never- 
theless, in both cases one finds that the escape probability scales 
like P ~ 1 - i Tj for t|j A Tj <k 1, but the scaling for intermediate 
cases can be rather different. 

Below we now will give the solution for N' including feed- 
back and Hi continuum absorption in the no scattering approxi- 
mation. The equations are applicable to both hydrogen and helium 
recombination, however during hydrogen recombination there is no 
continuum opacity that can affect the evolution of the photon dis- 
tribution between the H i resonances significantly. 



5.2 Feedback between isolated resonances 

In the standard quasi-stationary approximation for an isolated res- 
onance one has to set 0|, = 1 and = 1 in Eq. ( I4bb . This then 
yields FJ,'' = 1 - e"'"abs*''~s~'. Furthermore, one should assume that 
<t,(v,z',z) = T'^lx'y, -xil withal, = £'>f'(v')dv' and t'^ = p'^T'^, 
where Tj is the standard Sobolev optical depth. Inserting this into 



- q 

A'' ,. However, the optically thin limit is reached (slightly) (see lChluba & Sunvaevll2009d) 



Eq. {175) for -» oo one then has P^-^ jz) = P'/z) = [1 - e"^j]/r;, 
(e.g. see also lChluba & Sunvaevir2009dl) . 

In this approximation the spectral distortion on the red side of 
the resonance i is given by 



ANL 



AN'_ 



A^;„(z)[l 



(23) 



which is constant with frequency. Note that for the quadrupole and 
intercombination lines the factor 1 - e^'^J < 1 is important. 

For Eq. l l23t it was assume that on the distant blue side (y » 
y, ) of the resonance / the spectrum is given by the CMB blackbody 
and that the only distortion is created by the line itself. If we also 
allow for some distortion, AA'^,,, on the very distant blue side of the 
line then we have 



ANL 



AA'L(z)[l 



+ AN' e- 



(24) 



If one assumes that v, » v,_i then under quasi-stationary conditions 
the approximate behavior of the spectral distortion in the vicinity 
of the next lower-lying resonance y = ; - 1 is given by 

AWf = AA'^„,(z)[l - + AN'Az) e-<^'-'''\ (25) 

with the limiting cases 

„ . ( AN'Jz) + AN' (z) e-^i for v < y , 
ANf-i = \ (26) 
I AA'i(z) for V » Vj 

For this solution it was assumed that the cross-talk between the res- 
onances is negligible, meaning that in the region around resonance 
i the contribution of the opacity from i can be omitted. Note also 
that in our formulation AN'_{z} = AA'i(z). 



ANUz) 



^n[(z) - pi^ANJz) 



Pi 



(29) 



Inserting this into ANJz) = AA'en,(z) x [1 - P^^] and solving for 
ANJ{z) one then finds 



AA''(z) : 



1 - P\ 

' AA';(z) + 



pj 

'cm 



v' P' 
" " AN'_(z) 



^-pLPi 



AN[{z) - 



PiP, 



^ /-'em ^ A 



[aa'^(z)-aa':(z)1 



Pi Pi 

i-pLPi 



AN'Sz) 
AN[(z) 



(30a) 



(30b) 



(30c) 



Here P^'J^ is the effective Sobolev escape probability which can 
be directly used in the multi-level code and includes the effect 
of feedback. However, AN'_(z) still depends on AA'^,^(z), but with 
Eq. iSOaj it is now also possible to obtain an explicit expression 
for AA'em(z). which then also can be applied to compute AA'^^(z) 
knowing AA'!_+'(z). We find 



AA'em(z) = AA'£(z) + 



[a<(z)-AM+'(z)]- (31) 



For p*„, « 1 the second term can be neglected. Note that this 
limit can be reached for p*,^ « 1 and/or Pj « 1 . However, when 
Pj ~ 1 the latter term should be included. In particular even without 
feedback (AA'*;+'(z) = 0) one will have 



A<(z) {^^(Z) forpL/'*«l 

Ayv* (z) = — X I 

™ 1-pLP^, U^Nliz) forP*^l. 



(32) 



This shows that for < 1 one will find ANl(z) < AN^^{z). 
This fact is important in relation with the H i Lyman series and the 
He I n'Pi - I'So series during the pre-recombinational epochs of 
the considered species. 



5.3 Delay between emission and feedback 

As a next step we also take into account that the feedback of pho- 
tons from the resonance ; on j occurs at slightly different redshifts. 
This implies that for -Pf^'j at redshift z one has to evaluate AN'_ at 
z, = — [1 -I- z] - 1 > z. In addition one should take the change in 



the volume element into account so that 



A«/(c) 



AA'ifei) 
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>7 AJV/(;) 



A;i,'(c) 



, where Ara 



2v? 



is the deviation of the pho- 



ton occupation number from the one of the CMB blackbody. Then 
one finally finds 



PUz) 



p' 

d.S 



(z) 1 



1 



/''em ^ tl 



(33a) 
(33b) 



With this approximat ion also the feedback in th e H i Lyman-series 
was solved earlier fchluba & SunvaevI |2007| : ISwitzer & Hiratal 
l2008ah . Note that for the H i Lyman-series one can neglect the fac- 
tor 1 - e~^<i since tJj » 1 at all times during hydrogen recombina- 
tion. Also since t|j » 1 for all important H i Lyman-series transi- 
tions one will have An'_(z) An|.^j(z) x An'^^{z) and /"jjCz) ~ Piiz), 
so that pII^';.(z) ^ Pi(z)[l - A<(z,)/A«[(z)]. 

Because during the recombination of hydrogen the whole Ly- 
man series is extremely optically thick, all the photons released in 
the transition / will be reprocessed in the resonance j with vj < v, . 
In the case of the Hei «'Pi - I'So intercombination lines and 
He: ;!'D2 - I'Sq series (see Sect.|6]l this approximation is not al- 
ways justified, since the effective absorption optical depth for the 
Hei n^Pi - I'So transitions with n > 2 and at high redshift also 
for the Hei ?7'D2 - I'So series does not exceed unity at any red- 
shift (e.g. see Fig.[2]l. Therefore the feedback will not be restricted 
to i — » i - 1, but in some cases (a large) part of the distortion will 
also feedback like / — > i - 2. Also during the pre-recombinational 
epochs of helium and hydrogen one will have to use the full expres- 
sion Eq. l |33l l, since even some of the main resonances (e.g. the H i 
Lyman series) can become optically thin. 



5.4 P^^ J (z) and An* (z) for small Td 

In the pre-recombinational epochs of hydrogen and helium, and 
during the end of recombination one can have lines with rj « 1 
while Ts can in principle take all values. In particular for the Hi 
Lyman a line, because of the strong dependence of on redshift, 
one encounters the situation when tj <k 1 and ts » 1 . For « 1 
one can approximate Pi ~ 1 - | Td. Then one finds 



pisiz) 



. 1-4 



1 T 2 Ig 



(34a) 



and from Eq. i3l\ with Eq. i23l including the delay between the 
emission and feedback redshift one has 



=1 1-4 



An*(z) « Ti ^ 



1 i- 2 Ij 



A<(z)-pLA«-^'(zt+i)}. (34b) 



Here in addition the two limits for p*^Ts » 1 and /?*^rs ^ 1 exist, 
resulting in 

[1-|t* forr* « 1,p*„/s « 1 

Pas(z)«\ r (35a) 
[^[l-f| forr*«l,pL/s»l 



and 

A/ (z) « {Aniiz) - pLA//"'(z*+i)} 

""IeI^"^] f°>-r*«l,pLr*»l. 

(35b) 

In particular we can see that in the optically thin limit (t^ « 1 
and p*nTj « 1) the effective escape probability and the occupation 
number on the red side of the resonance both approach the values 
expected from the Sobolev approximation. 



5.5 Feedback between isolated resonances in the presence of 
neutral hydrogen 

While He i photons redshift from the resonance i to the next lower- 
lying resonance j = i - I there is also some absorption in the H i 
Lyman continuum. In order to take the additional absorption into 
account we assume that between the resonances the photons only 
feel the Hi continuum opacity. Once the photons, starting at reso- 
nance ( with An'^, X An'_ = An[.,^(z,) [l - e^'^j'^'^j, will have reached 

the resonance j at redshift z, one will have An\, x An'_(z,)e "^"J, 
where we approximate the continuum optical depth by 



ccri.sc(V;)A'"' V, - Vj 



H 



(36) 



More accurately one can use t':^j(z) = ■r!j|,s(vj, z,, z), where 
T^jj^(vy, z,, z) has to be computed using Eq. ilbl . This leads to 



□Sfc; 



1 - 



An'Jz,) ^ 
Aniiz) ' 



(37) 



A simil ar approximation was also given earlier bv lSwitzer & Hiratal 
l l2008al) . In the case of He I 2 ' P 1 - 1 ' So ^ He I 2^ Pi - 1 ' So feedback 
it works very well (cf. Fig.[6]and see Sect. l6.3t . 



5.6 Feedback on the H i Lyman continuum 

To include the feedback of Hei photons on the Hi Lyman con- 
tinuum, we follow a very simple procedure. We assume that ev- 
ery resonance produces a distortion to the photon occupation num- 
ber of A;iJ„(z') = An'_(z') at v,. Due to redshifting this distortion 
then moves towards lower frequencies, so that at redshift z and 
frequency v one will have a distortion of An'^(z) ~ An'_{zi) with 
z, = 7 [1 + z] - 1 > z. If we now take into account that on their 
way some of these photons are also absorbed in the Hi Lyman 
continuum then we will find An'^(z) ~ A«'_(z,) e"^*s*''-'-', where 
'^IbS^'Zi'Z) has to be computed using Eq. i7b\ . Then the correc- 
tion to the ionization rate of the H i Is state caused by the distortion 
from the resonance / is approximately given by 



AR\^ « 47r tri,,(v) —An'_{z,) e-^-bJ"- :) dv. 



(38) 



This integral can be computed for every line during the run of the 
multi-level recombination code. Note that whenever the photons 
from resonance ; are passing through a resonance j < i then the 
distortion at v < is suppressed by e^'^J. For the Hei n'Pi - I'So 
series and at high redshift also for the Hei n'Dj - I'So series this 
strongly increases the distortion at frequencies below the line, a fact 
that has to be taken into account in the computations. 
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Figure 7. Escape probability in the He i 1 ' So and H i Lyman continuum. 
The curves were computed using the approximation for the escape proba- 
bility given in Appendix lB2l 



5.7 Feedback by photons escaping the H i Lyman continuum 

Even though the direct escape of photons in the H i Lyman contin- 
uum is not important during the epoch of hydrogen recombination, 
in the pre-recombinational epoch of hydrogen (z > 1800) the es- 
cape probabihty can be close to unity (cf. Fig.[7]l. Due to the feed- 
back of helium photons on the H i Lyman continuum one can there- 
fore produce some amount of photons that directly escape below 
the H I Lyman continuum frequency. These photons then feedback 
in the forrest of H i Lyman series lines, where for the uppermost 
transitions one could in principle model this process as a contin- 
uation of the H i con tinuum cross section (e.g. see Appendix II in 
lBeigmanetal]|l968l) . Here we approximately accounted for this 
process by simply adding the small distortion from the H i Lyman 
continuum on the blue side of the uppermost feedback level that we 
included. The distortion from the H i Lyman continuum was com- 
puted using the net rate as given by Eq. dlSK and assuming that 
the photons were emitted at the Hi Lyman continuum threshold 
frequency. 

In addition, some of the high frequency photons emitted in the 
pre-recombinational epoch of He i induced by the feedback of He ii 
photons on neutral helium can reach below the Hi Lyman contin- 
uum frequency. We also included these photons in the computations 
by feeding them back on the blue side of the uppermost feedback 
level in hydrogen. 

For the processes in connection with He ii — » He i feedback we 
followed a similar approach. However, in this case the distribution 
of photons from HEin — > Hed recombination can be computed 
before solving the problem (see Sect.[8]for details). 
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6 CORRECTIONS TO THE IONIZATION HISTORY 
DURING HELIUM RECOMBINATION 

Below we will now look at the feedback of photons between the 
different resonances of helium and their effect on the ionization 
history. Our main results are given in Fig.|9]and[T0] At low redshift 
our final correction seems to be smaller than the one presented in 



Figure 8. Corrections to the free electron fraction during He ii — ► He i re- 
combination: effect of additional Hei n'Pi — I'Sq transitions (upper panel), 
Hei n'Di - I'So electric quadrupole transitions (middle panel), and the 
effect of the hydrogen continuum opacity on the Hei n'Pi — I'Sq series. 
We included 10 shells for hydrogen and 10 shells for helium in our com- 
putatio ns, and compared the resul ting electron fraction with our reference 
model iRubiiio-Martm et alj20O8l) . 
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Table 1. First few main transitions to the ground-state of helium. The res- 
onances have been ordered according to their transition frequency. In this 
way one can understand the sequence of feedback between the different 
lines. For this purpose we also give the A;/z ~ Ay/y that is necessary to 
reach the next lower-lying resonance. For comparison due to the thermal 
motion of the atoms one has Ayo/y ~ 1.7 X 10"^ [^m] 



Initial Level 


Type 


Ails J 


>'ils LnZJ 


1 A - /-I 






El (TS) 


177.58 


5.0691 X 10'5 






Tip, 


CI /CCA 


1 TOSQ \y 1 a9 
i . lyoy X iU 


^ 1 \/ 1 nl5 

J. 1 jUj X iU 


i.Z X lU 


-2 


33P, 


El (TS) 


56.1 


5.5631 X 10'5 


7.8 X 10- 


-2 


3'D2 


E2 (SS) 


1.298 X 10^ 


5.5793 X 10'^ 


2.9 X 10" 


-3 


3'Pi 


El (SS) 


5.6634 X 10** 


5.5824 X 10'^ 


5.6 X 10- 


-4 


43pi 


El (TS) 


23.7 


5.7325 X 10'5 


2.6 X 10" 


-2 


4ID2 


E2 (SS) 


748.5 


5.7394 X 10'^ 


1.2 X 10- 


-3 


4'Pi 


El (SS) 


2.4356 X 10** 


5.7408 X 10'^ 


2.4 X 10- 


-4 


53P, 


El (TS) 


12.1 


5.8100 X 10'5 


1.2 X 10- 


-2 


5'Do 


E2 (SS) 


431.4 


5.8135 X 10'^ 


6.0 X 10- 


-4 


5'P, 


El (SS) 


1.2582 X 10** 


5.8143 X 10'^ 


1.4 X 10- 


-4 



ISwitzer & Hiratal l l2008ij) . However, it is clear that the difference 
will not be very important for the analysis of future CMB data. 
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Figure 9. Corrections to the free electron fraction during He n ^ He i re- 
combination: the effect of feedback among the resonances of neutral he- 
lium. We show the results for two sets of computations, one that does 
not include the effect of photon absorption between the lines (labeled 
'no LL-HI opacity') and those that do ('with LL-HI opacity'). We in- 
cluded 10 shells for hydrogen and 10 shells for helium in our computa- 
tions, and compared the resu lting electron fraction with our reference model 
)Rubiiio-Martin et al.ll2008l) . 



6.1 Initial reflnements of the recombination code 

In this Section we mainly want to discuss the effect of feed- 
back of helium photons on the cosmological recombination his- 
tory. However, for this we also need to take into account additional 
processes, which we neglected until now. First we also include 
He I n' D2 - 1 ' So electric quadrupole transitions and He i n^Vy - 1 ' So 
intercombination lines with n > 2, which all have transition rates 
and escape probabilities that should lead to net rate which are com- 
parable to those of the main transitions controllin g the dynamics 
of helium recombination jSwitzer & Hiratal2008ah . In Appendix IaI 
we give some more details on the rates that we used. In Fig. [8] we 
show the effect of these transitions on the recombination history. 
The corre ctions that we find seem to be in good agreement with the 
results of lSwitzer&Hiratj ( l2008ah for these processes, however 
the final corrections are negligible. 

Then we also refine our modeling of the escape of photons 
from the Hei «'Pi - I'So series with n > 2. Like in the case 
of the He: 2'Pi - I'So series and the Hei l^Fi - I'So the pres- 
ence of neutral hydrogen increases the escape probabilities of these 
lines towards the end of helium recombination. We take this pro- 
cess into account using th e ID-integral approximation given by 
iRubino-Martm et al. (see Eq. (B.3) in their paper) . We show 

the effect of this proces s on the ion ization history in Fig. [8] Again 
in good agreement withj Switzer & H irata (2008a) we find a negli- 
gible correction. We do not include this speed-up for the escape of 
photons in the Hei ra'D2-l'So and Hei ra^Pi-l'So series forn > 2, 
since the effect of these lines on the dynamics of helium recombi- 
nation is already very small. Additional refinements of the escape 
probabilities will not change this very much, sin ce the speed-up in 
the mai n resonances will dominate. The results of lSwitzer & Hiratal 
( l2008ij) also support the correctness of this statement. 



6.2 Total amount of primary He i photons which are 
available for feedback 

As shown b\ l^lRubino-Martfn et al ](|200i), about 90% of all helium 
atoms became neutral via the Hei 2'Pi - I'Sq and Hei 2'Pi - 
I'So dipole transitions to the ground state. Direct transition to the 
ground state from initial levels with principle quantum number n > 
2 allowed about 2% of the helium atoms to recombine, while the 
remaining ~ 8% of helium atoms became neutral via the He i 2' Sq- 
1'Sq two-photon decay channel. This implies that in total about 
1 .08 7 per helium atom were emitted in these transition, where ~ 
1 7 per helium atom was released above the H i Lyman a resonance. 

It is clear that the photons emitted in He i A'-'^ + ' Ly - 1 ' So tran- 
sitions all have the possibility to feedback on the H i Lyman contin- 
uum, while a large part of the Hei 2' So - I'So two-photon contin- 
uum will feed back only on the H i Lyman series (see Sect. 17.21 ). 

6.3 Hei2'Pi - I'So ^ HEi23Pi - 1' So feedback 

In Table[T]we give the main resonances, that are important for the 
feedback problem. Here we first consider the feedback of photons 
from the Hei 2'Pi - I'So line on the Hei 2''Pi - I'So resonance. 
These resonances are separated by Av/y ~ 1% or ~ 700 Doppler 
width, so that one can certainly treat them as two isolated reso- 
nances. For this rather large distance one should in principle com- 
pute the H I continuum opacity between the lines with the full inte- 
gral Eq. l |7bt rather than the approximation Eq. i36\ . However, for 
the accuracy required here we still use the simple approximation. 



It turns out that also with the refinements to our computations introduced 
in this paper the quoted numbers do not change very much. 
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since otherwise for the current version of our code the computa- 
tions are slowed down significantly. 

In Fig.|9]we show the results of our numerical computations 
for this problem. If we neglect the absorption of photons in the 
Hi continuum while they travel from the Hei 2' Pi - I'So line to- 
wards the Hei - I'So resonance we obtain the thick dashed 
curve. This leads to a significant correction during helium recom- 
bination, which is largest at z ~ 2020 with ANJN^ ~ -1-0.36%. 
However, when we include the effect of the continuum opacity (thin 
dashed line), a large part of the correction disappears, leaving only 
ANJN^ ~ +0.1% at z ~ 2300. Looking at Fig.[3](sohd line) this 
behavior is not surprising, since at z ~ 2370 the Hi Lyman con- 
tinuum becomes optically thick between the two resonances. This 
stops Hei 2'Pi - I'So photons from reaching the Hei 2'Pi - I'So 
resonance. This effect also can be seen in Fig.|6] where at z <; 2600 
the reduction in the escape probability from the Hei 2'' Pi - I'So 
line is only caused by the feedback of photons. At z ~ 2600 the 
efi^ect of the hydrogen continuum opacity starts to set in and until 
z ~ 2200 none of the photons from the He: 2'Pi - I'So ever reach 
the intercombination line. 
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Figure 10. Ionization history during He n — > He i recombination. The re- 
sults were obtained for 10 shells for hydrogen and 10 shells for helium 
including the feedback between al l 10 shells of helium. We also changed 
our cosmology to match the one o f ISwitzer & Hiratal j2008bl) . 



6.4 Feedback among higher level transitions 

In this Section we now include the effect of feedback from the 
higher level transitions. Here the interesting aspect is that for each 
shell n > 2 one has a feedback sequence Hei ;i'Pi - I'Sq — > 
He: n'D2 - I'So — > Hei Pi - I'So where the separation between 
the lines is rather small. For example, the He: 3 'Pi - I'So pho- 
tons will feedback on the He: 3'D2 - I'So transition after Az/z ~ 
5.6 X lO"'*, since the quadrupole line at z ~ 2200 is only about 30 
Doppler width below the Hei 3'Pi - I'So resonance. Similarly the 
He: 3^^Pi - I'So resonance is only 170 Doppler width away from 
the He I 3 ' D2 - 1 ' So quadrupole transition (see Table[T]for more ex- 
amples). In particular for the He: n'Pi - I'So — * Hei «'D2 - 1 'So 
feedback it is questionable if it is really possible to treat the reso- 
nances completely independently, since the line broadening due to 
resonance and (to a smaller extent) electron scattering will likely be 
larger than the separation of these l ines. For example in the ca se of 
hydrogen our computations show l lChluba & Sunvaevll2009al) that 
the line broadening typically is of the order of a few x 10^ Doppler 
width in the case of Lyman a. For the Hei «'Pi - I'So sequence 
it will be a bit smaller due to the fact that there is less helium than 
hydrogen, but still it should exceed few x 10 Doppler width. How- 
ever, we will still use our simple approximation, as the correction 
again turns out to be rather small. 

As mentioned above, since the He: 3^^Pi - I'So line is never 
really optically thick (cf. Fig. |2j, most of the photons from the 
He: 3'D2 - I'So quadrupole line will pass through this resonance 
and therefore could also feedback on the He: 2' Pi - I'So line. How- 
ever it takes about Az/z ~ 8% to travel from the Hei 3^Pi - I'So 
to the Hei 2'P, - I'So line. Looking at Fig. [3] we can see that at 
z ^ 2600 - 2700 photons emitted by transitions with /; = 3 will 
never reach the He: 2'Pi - I'So resonance, because they will be re- 
absorbed in the H i continuum before. Similar comments also apply 
for the higher Hei n^Pj - I'So series and at early times even for 
the photons from the Hei n'Di - I'So quadrupole transitions. For 
example, because of the small fraction of neutral hydrogen atoms 
present during Hei recombination photons emitted in transitions 
n = 4 to )2 = I will no longer reach the He: 3 'Pi - I'So resonance 
at z < 2300, and at z < 1900 - 2000 Hei 3'Pi - I'So photons 
will not be able to feedback on the He: 3'D2 - I'So resonance 



(cf. Fig.[3j. For our computations we also took these aspects of the 
problem into account. 

In Figure |9] we also show the results of our computations in- 
cluding the feedback for the higher transition, iif = 3 means that we 
took the feedback between the main resonances with « < 3 into ac- 
count, starting the sequence with the He i ;i ^ P 1 - 1 ' So line . Again we 
considered the cases with and without including the hydrogen con- 
tinuum opacity between the subsequent resonances. The reabsorp- 
tion of photons between the lines leads to a very large suppression 
of the feedback correction, which for ;if = 10 reaches AN^/Ne ~ 
0.17% at z ~ 2300 instead of ANJN^ ~ 0.46% at z ~ 2045. We 
also ran cases including more than 10 shells in the computation, and 
found that the correction already converges for rif ~ 10. However, 
given that the final result is rather small we did not investigate this 
in more detail. One can also see from Fig.|9]that the largest single 
contribution to the feedback is coming from the second shell, i.e. 
the feedback sequence Hei 2'P, - I'So He: 2^Pi - I'So. 

In Fig. |10 | we also give the direct comparison with the results 
presented in Iswitzer & Hiratal ( l2008bl) . Our feedback correction 
seems to be slightly smaller at z g 2400, when including the effect 
of photon absorption in the Hi Lyman continuum between subse- 
quent resonances. Su rprisingly, our result seem s to be very com- 
parable to the one of ISwitzer & Hiratj JioOSbl) , when we do not 
include this additional continuum absorption. However, when in- 
cluding the H I continuum opacity between sub-sequent resonances, 
our result is slightly smaller. Nevertheless one should mention that 
these kind of changes in the helium recombination history will not 
be very important for the analysis of future CMB data. 



6.5 Feedback of He i photons in the H i Lyman continuum 

We also computed the recombination history including the feed- 
back of He I photons on the H i Lyman continuum, using the ap- 
proximation discussed in Sect. l5.6l As we already mentioned above 
all Hei photons are feeding back on hydrogen during its pre- 
recombinational epoch. Therefore one does not expect any large ef- 
fect in the ionization history. We found that the correction does not 
exceed AN^ /N^ ~ 1 0"* - 1 0"^ so that for computations of the CMB 
power spectra this process is negligible. However, as we show in 
Sect.|7] this process does lead to some traces in the cosmological 
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Figure 11. Effect of Lyman series feedback on the ionization history during 
H I recombination. The results were obtained including 20 shells for hydro- 
gen and 20 shells for helium and accounting for the feedback between the 
first 10 shells of hydrogen. 



recombination radiation, which should still be present as distortion 
of the CMB today. 

We also checked whether the feedback of the helium induced 
pre-recombinational features from the higher H i Lyman series will 
have an effect on the dynamics of h ydrogen rec ombination. Like in 
the normal recombinational epoch Jchluba & S unyaev 200^) these 
additional photons will feed back on the next lower-lying Lyman 
series transition, however, due to the fact that at high redshift the 
Lyman series is not completely optically thick, the feedback is no 
longer restricted to Ly;i — > Ly(« - 1), but can go beyond that. We 
refined our computation of the H i Lyman series feedback problem 
including this process, but found no significant correction during 
the pre-recombinational epoch. 

Furthermore, we also included the small feedback correc- 
tion due to photons that are emitted in the Hi Lyman continuum 
during the pre-recombinational epoch (see Sect. 15.71 for more de- 
tails). Again the correction to the ionization history was negligi- 
ble. However, as we explain below (Sect. |7j this process leaves 
some interesting traces in the recombinational radiation spectrum. 
Our final result for the effect of feedback during hydrogen recom- 
bination is shown in Fig. [TT] Note that the result is a bit bigger 
(by about AN^/AN^ ~ .04% at z ~ 1050) than presented in 
IChluba & SunvaevI ilOO'k because we slightly improved the accu- 
racy of our numerical treatment for the normal Hi Lyman series 
feedback. 



7 PRE-RECOMBINATIONAL EMISSION FROM 

HYDROGEN DUE TO THE FEEDBACK OF PRIMARY 
HE I PHOTONS 

In this Section we present the results for the cosmological recom- 
bination spectrum from hydrogen and He i, when including the full 
feedback of primary He i photons from the «'Pi - I'So, """Pi - I'So 
and «'D2 - I'So series with n < 10. Figure[T2l shows the results of 
our computations. Starting with the helium bound-bound contribu- 



bouiid-bound HI recombination spectrum 

normal recombination spectrum 
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normal recombination spectrum 
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Figure 12. Contributions to the cosmological recombination spectrum. We 
show the H i bound-bound spectrum (upper panel), the H i free-bound spec- 
trum (middle panel), and the He: bound-bound spectrum (lower panel). 
We included 20 shells of hydrogen and 20 shells of helium in our com- 
putations and allowed for helium feedback from the first ten shells. The 
curves marked with 'RMCS 2008 ' where computed using our reference 
model iRubino-Marti'n et alj|2008l ). Those marked as 'normal recombina- 
tion spectnim' do not include the effects of the hydrogen continuum opac- 
ity on the dynamics of helium recombination. Note that for the He i bound- 
bound spectrum we plotted the modulus of the distortion. There are two 
negative features at v ~ 140 GHz and v ~ 270 GHz. 
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tion one can see tha t not very much has changed with respect to our 
reference spectrum jRubino-Martm et alJlOOSl) when including the 
He I feedback process. The most important difference is that there is 
no high frequency spectral distortion due to direct transitions from 
levels with > 2 to the ground state. These photons have all died 
in the H i Lyman continuum during the pre-recombinational epoch 
of H I, and also for our reference spectrum we neglected them here. 

If we look at the bound-bound contribution from hydrogen 
(Fig.[T2l upper panel), then we can see that the feedback of helium 
photons leads to a pre-recombinational H i feat ure at v ~ 1300 GHz. 
Part o f this feature was already given by iRubino-Martm et aD 
( l2008l) . however, due to the additional reabsorption of high fre- 
quency photons from neutral helium, this signature of helium feed- 
back increased. As we will see below the total amount of pho- 
tons appearing in that feature practically doubled. Note, that also 
at lower frequencies the helium feedback leads to some notable 
changes in the bound-bound emission from H i. 

In Fig. [T2] we also present the changes in the Hi free- 
bound contribution, which were not presented elsewhere so far. 
One can see that here no strong additional helium induced pre- 
recombinational features are visible. This is because the Hi free- 
bound continua are rather broad at all times (typical Av/v, ~ 
kTy/hvj ~ 0.17 X ^) and therefore overlap strongly, so that 
He I feedback only leads to an increase in the overall amplitude of 
the free-bound continuum. We also checked the emission in the 2s- 
Is two-photon continuum, but the changes were so small that we 
simply neglected this contribution here. 



7.1 Counting the number of additional photons 

We can now look at the number of additional photons that are cre- 
ated due to the feedback on hydrogen. It is clear that every helium 
photon that is absorbed in the Hi Lyman continuum will at least 
be replaced by one photon in the free-bound spectrum of hydro- 
gen. Since so far we only took the feedback of photons emitted in 
theHEi n'Pi - I'So, Hei n^Pi - I'So and He: n'Dj - I'S,, series 
with n < 10 into account, and because these amount to about 0.9 
photons per helium atom available for the feedback process, one 
expects that the free-bound emission from hydrogen increases by 
about 7%. Simply computing the total amount of photons emitted 
in the free-bound continuum of hydrogen, we indeed find that in 
total ~ 7.1% more photons are created in the Hi free-bound con- 
tinuum (Table [2J. About half of this number is already appearing 
for our reference model, in which only part of the feedback correc- 
tion is included. 

We then also looked at the number of photons emitted in the 
H 1 2s- 1 s two-photon continuum, finding that this number in not af- 
fected much due to helium feedback (cf. Table[2j. This is because 
at early times the Lyman a line and the other Lyman series transi- 
tions are more optically thin and therefore mu ch more efficient than 
the 2s - Is two-photon channel. As shown in iRubino-Martm et al.l 
( l2006h the 2s- Is two-photon channel only starts to become active 
at redshift z ;£ 1600 - 1700 where practically no additional helium 
photons are produced. 

We also computed the total number of additional photon emit- 
ted in the Lyman series and found about 7% per hydrogen atom. 
Together with the contribution from the 2s- Is two-photon channel 
this yields ~ 7.1% per hydrogen atom as it should be, since all 
electrons that entered the hydrogen atom via the continuum, should 
reach the ground state via the Lyman series or alternatively the 2s- 
Is two-photon channel. From Table [2] one can also see that most 



Table 2. Total number (density) of photons, Ny, produced by different tran- 
sitions in hydrogen. Here A'h and Whc denote the total number (density) of 
hydrogen and helium nuclei. Contributions marked with 'rec' are from the 
normal recombinational epoch. Those mai'ked with ' pre-rec F are the ad- 
ditional photons computed for our reference model iRubiiio-Martm et alj 
2008). while those marked with 'pre-rec IF give the additional photons 
when accounting for the full He i feedback. We also give the numbers for the 
whole bound-bound and free-bound spectrum, as well as the total difference 
(bound-bound + free-bound) in the number of photons. 



Line 






H I hya (rec) 


42.5% 


— 


H I hya (pre-rec I) 


3.3% 


42% 


H 1 hya (pre-rec II) 


3.1% 


39% 


H I hya (all) 


48.9% 




H I Ly series (rec) 


42.7% 


— 


H I Ly series (pre-rec I) 


3.6% 


46% 


H I Ly series (pre-rec II) 


3.4% 


43% 


H I Ly series (all) 


49.7% 




Hi 2s- Is (rec) 


2 X 57.2% 


— 


Hi 2s- Is (pre-rec I) 


~ 


— 


Hi 2s- Is (pre-rec II) 


2 X 0. 1 % 


1% 


Hi2s-ls(aU) 


2 X 57.3% 




H I free-bound (rec) 


1 nCict, 




H I free-bound (pre-rec I) 


j.O /c 


to /o 


H I free-bound (pre-rec II) 


3.5% 


44% 


H I free-bound (all) 


107.1% 




Hi bound-bound (dipole, low-y, rec) 


1.94 




H I bound-bound (dipole, low-y, pre-rec I) 


5.3% 


67% 


H I bound-bound (dipole, low-y, pre-rec II) 


5.1% 


65% 


H I bound-bound (dipole, low-v, all) 


2.04 




H I bound-bound (dipole, rec) 


2.37 




H I bound-bound (dipole, pre-rec I) 


8.9% 


113% 


H I bound-bound (dipole, pre-rec II) 


8.5% 


108% 


H I bound-bound (dipole, all) 


2.54 




total difference (pre-rec I) 


12.7% 


1.6 


total difference (pre-rec II) 


12.2% 


1.6 


total difference (pre-rec I+II) 


24.9% 


3.2 



(i.e. ~ 6.4%) of these additional high frequency photons are pro- 
duced in the Lyman a line, but also several additional photons are 
created in the higher Hi Lyman-series. 



7.1.1 Total number of photons related to He i 

We finally also computed the total number of additional photons 
which are related to the presence of neutral helium in our Universe 
(see Table|2ll. Per Hei atom we found about 3.2 additional photons 
in the H i recombination radiation. Since in total about ~ 0.9 ioniz- 
ing quanta per helium atom were absorbed in the feedback process, 
and then eventually replaced by a Hi Lyman a resonance or 2s- 
Is two-photon continuum photon (see Sect. I7.lt . about 2.3 extra 
photons per helium atom appeared in the H i recombination spec- 
trum. These photons would be absent without the inclusion of the 
He i feedback process. Furthermore, of these 2.3 extra photon per 
helium atom about ~ 0.9 y were produced in the free-bound contin- 
uum and the remaining ~ 1.4 y were emitted in transitions between 
excited states with principle quantum number n > 2. Furthermore, 
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one can say that per He i feedback photon about 2.3/0.9 ~ 2.6 extra 
photons were released by Hi. 

In our reference model per He i atom we find at total of ~ 1 .85 
photon^in the normal He i bound-bound spectrum, ~ 0.16 photons 
in the He i 2' So - 1 ' So two-photon decay spectrum, and ~ 1 photon 
in the He i free-bound continuum. This implies a total of ~ 3 y per 
He I atom. Including the full feedback caused by the photons from 
the Hei fi'Pi - I'So, Hei fi^p, - iiSo and Hei (i'Dj - I'So series 
with n < 10, we removed a total of ~ 0.9 high frequency photons 
per helium atom, but add ~ 3.2 new photons to the Hi recombina- 
tion spectrum. Therefore in total the number of photons related to 
the presence of He i increased to a total of 6.2 y per helium atom, or 
by a factor of ~ 1.7 relative to the case with no interaction between 
helium and hydrogen for which a total of 3.7 y per helium nucleus 
are emitted (see Table[3}. This is a large increase in the total number 
of photons that are related to the presence of He i in the Universe. 
One should also mention, that the absolute number should increase 
when including more shells for both the hydrogen and neutral he- 
lium atom into the computation. Also, the feedback from higher 
level transitions and the photons emitted in the He: 2' So - I'So 
two-photon continuum will slightly increase this number (see Sec- 
tions below). We leave a detailed computation for a future paper. 

7.2 Additional corrections to the He i feedback problem 

7.2.1 Feedback from higher level transitions 

Here we restricted our computation to the feedback of photons from 
the Hei n'Pi - I'So, Hei n^Vi - I'So and He: w'Dt - I'So series 
with n < 10. We computed the additional number of photons emit- 
ted in these series from n > 10 and found about 0.02 photons per 
helium atom. Taking the feedback of these into account one there- 
fore expects 8% x 2% x 2.6 ~ 0.4% extra photons per hydrogen 
atom (or ~ 0.05 y per helium atom). As we discuss below also the 
feedback from photons in the He: 2' So - 1 'So two-photon contin- 
uum should lead to another small increase in the total number of 
photons that are related to the presence of helium in the Universe. 

7.2.2 Feedback of photon from the He i 2' So - 1 ' So two-photon 
continuum 

In our computation so far we did not include the correction due 
to the feedback of photons emitted in the Hei 2' So - I'So two- 
photon continuum. It is clear that ~ 1 /2 of these photons are never 
able to feedback on hydrogen, since they are directly emitted below 
the Hi Lyman a line with transition energy 10.2 eV. Therefore, as 
mentioned in the Sect. |6.2| in principle per helium atom there are in 
addition ~ 0.08 photons available for the feedback on hydrogen. Of 
this number about 0.044 photon^ per helium atom will feedback 
in the hydrogen Lyman continuum, while the remaining photons 
could directly affect the H i Lyman series. Again this feedback will 
not introduce any important change in the ionization history, but 
will only lead to the production of additional H i photons. Assum- 
ing that all these photons are absorbed one should therefore find 
another 8%x8%x2.6~ 1.7% more photons per hydrogen atom, or 
~ 0.21 additional photons per helium atom. The final number is ex- 
pected to be a bit smaller, because not all photons will be absorbed 

' In this number we neglected the high frequency resonance photons that 
will eventually be reprocessed by H i in the feedback problem. 

We estimated this number by integrating the He i 2' So- 1 ' So two-photon 
profile over the appropriate range of frequencies. 



and since the feedback in the Hi Lyman series will not produced 
as many secondary photons, but we leave a detailed calculation for 
some future work. 



8 FEEDBACK OF PHOTON FROM HE III ^ HE II 
RECOMBINATION 

It is clear that also the feedback of high frequency photons from 
He III — > He ii recombination will lead to some additional changes 
to the CMB distortions from cosmological recombination. Like in 
the case of He i feedback on hydrogen, there is no significant addi- 
tional corrections to the ionization history caused by this process, 
since He ii feedback is again occurring in the pre-recombinational 
epoch, but this time of He i. However, several secondary photons 
are produced, increasing the total contribution of helium related 
photons to the cosmological recombination spectrum. A rigorous 
computation of this problem is beyond the scope of this paper, 
however, there are several interesting aspects, that we would like 
to mention, before demonstrating the principle importance of this 
feedback process for the total contribution of helium-related pho- 
tons to the cosmological recombination spectrum (see Sect. l8.8t . 

8.1 Total amount of high frequency He n related photons that 
are initially available for feedback 

As mentioned in the introduction, all He ii photons emitted at ener- 
gies abovj*^ ~ 20.96 eV (corresponding to the transition energy of 
the Hei 2^P| - I'So resonance) in principle can feedback on Hei, 
while those photons emitted at energies above ~ 10.2 eV (i.e. the 
Hi Lyman a transition energy) may affect hydrogen. These state- 
ments are rather simple, but the details turn out to be more involved. 
Below we provide some estimates for the total number of He ii pho- 
tons that are available in the feedback problem. 

8. 1. 1 Initial amount of high frequency He ii photons that are 
available for the feedback on He i 

According to our computations, for Heii recombination about 45% 
of Heii electrons went through the Heii 2s- Is two-photon contin- 
uum, while the remaining ~ 55% recombined via the He n Lyman 
a channel (with a very small addition due to the higher Heii Ly- 
man series). Furthermore, about 0.2 y photonj^per helium nucleus 
were emitted in the He ii Balmer continuum, and in total ~ 4.69 y 
were produced per Heii atom. These numbers were obtained in- 
cluding 100 shells for the Heii atom. 

All the Heii Lyman a photons (with energy ~ 40.8 eV) are 
able to directly ionize Hei atoms. After significant redshifting, 
these photons in principle could also directly ionize hydrogen 
atoms, but it turns out that all Heii Lyman a quanta will already 
disappear in the He i I'So continuum before that (see Sect. |8.6| l. To 
simplify matters, for our more detailed estimates of the effect on 
the cosmological recombination spectrum (Sects. [8j9j we will only 
take these photons into account. 

It is clear that also the Heii 2s- Is two-photon continuum will 

' ' In principle He ii photons could also feedback on the He i 2 ' So — 1 ' So 
two-photon continuum, but this transition is never opticaffy thick during 
recombination so that we neglected it here. 

T his number was estimated f rom the result for the 100 sheif hydrogen 
atom iChluba & Sunvae\j|2006ah . 
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Figure 13. Effective absorption optical depth, = pJjTj, in the Hi Lyman 
series. An interpretation and definition of rjj is given in Sect. 13. II 

lead to some interesting feedback. Abou{^ 37% of the photons 
emitted in the Heii 2s- Is two-photon continuum have energies 
above the Hei I'Sq continuum. Therefore, in total one should have 
~ 0.33 7 per helium nucleus with energy larger than the He i ion- 
ization potential E ~ 24.59 eV. 

In addition, those photons emitted in the Heh 2s- Is two- 
photon continuum with energies 20.96 eV < Ay < 24.59 eV (i.e. 
about 0. 1 7 per helium nucleus) may be able to directly feedback 
to the Hei )i'P, - I'Sq, Hei n^'Pi - I'So and Hei ?!'D2 - I'So se- 
ries, depending on the optical depth in these resonances at different 
stages during recombination. However, since the Hei 2'Pi - I'So, 
Hei 2'Pi - I'So and Hei 2'D2 - I'So resonances are not optically 
thick at z > 3000 (see Fig.|2j, most of these Heh photons will sur- 
vive and instead be available for feedback onto hydrogen at lower 
redshift. Also it is clear that the Heh Balmer continuum photons 
will not feedback on He i, but could affect the recombination of hy- 
drogen (see next Section). 

&.1.2 Initial amount of high frequency Heii photons that are 
available for the feedback on H i 

Regarding the possibility of direct ionization of hydrogen atoms by 
high frequency He n photons one will have those photons from the 
Heh 2s- Is two-photon continuum emitted in the frequency range 
13.6 eV < hv < 20.96 eV, amounting to about 0.21 7 per helium 
nucleus. Furthermore, those Heh 2s-1s photons emitted at energies 
10.2 eV < hv < 13.6 eV (~ 0.09 7 per helium nucleus) will be 
able to feedback on the Hi Lyman series. And finally, all those 
photons emitted in the Heh Balmer continuum in principle will 
also be able to feedback on the Hi Lyman series. This implies a 
total of ~ 0.5 potential feedback photons per helium nucleus (or 
0.04 photons per hydrogen nucleus) that may directly feedback on 
hydrogen. Even though these photons should lead to secondary H i 
Lyman series photons in the pre-recombinational epoch which then 



We estimated the numbers given below by integrating the Heh 2s- Is 
two-photon profile over the appropriate range of frequencies. 



again feedback on hydrogen itself, here we will not consider this 
part of the feedback process any further. 

Also the numbers given above should be considered as esti- 
mates since at z > 2500 none of the Hi Lyman series transitions 
is optically thick (cf. Fig. I13t , implying that the feedback process 
will not re-process all these photons. Similarly, some photons may 
escape the feedback problem in neutral helium and therefore add 
to the total amount of primary feedback photons available for hy- 
drogen. In addition, practically every ionizing He ii photon that is 
re-processed by He ii will be replaced by some secondary high fre- 
quency He I photon (in addition to some secondary low-frequency 
photons that do not affect the recombination problem anymore). 
This leads to a double-reprocessing of He n photons (see Sect. |8.2| l. 
A detailed analysis is beyond the scope of this paper. 

8.2 Double-reprocessing of He n photons 

The probably most interesting aspect of the Heu feedback prob- 
lem is the double-reprocessing of Heii photons, first by neutral he- 
lium and then later by hydrogen. This is because in addition to the 
primary photons initially released by He ii there will be secondary 
photons produced by both He i and H i. Those from hydrogen will 
only feedback on hydrogen itself but likely will not produce many 
more (low frequency) photons. However, the feedback of secondary 
Hei photons on hydrogen is expected to be more important. The to- 
tal number of these secondary He i photons will be very close to the 
total number of He ii photons that were re-processed. This can be 
understood because in the pre-recombinational epoch of He i the 
Hei 2' So - I'So two-photon channel will not be able to compete 
with the main resonances connecting to the ground state. Therefore 
it is clear that practically every He n feedback photon, independent 
of the way it was re-processed, will be replaced by a high frequency 
He i photon, which then again could be re-processed by hydrogen. 
Similarly, every high frequency photon that is re-processed by hy- 
drogen will eventually be re-appearing as a H i Lyman a photon. 

Using the numbers given in the previous section one there- 
fore expects a maximum of ~ 0.98 primary feedback photons per 
helium nucleus, and ~ 0. 12 primary plus secondary feedback pho- 
tons per hydrogen atom that are related to He n. 

In every feedback there is the potential possibility to produce 
several additional photons at low frequencies because electrons 
that reach highly excited states (via excitations from the ground 
state or direct recombination after feedback) can cascade towards 
lower states via many intermediate levels, emitting low frequency 
photons on their way. Therefore this process can lead to an im- 
portant increase of the helium-related photons in the cosmological 
recombination spectrum. 

8.3 Direct recombinations to the ground state of He i and H i 
and the efficiency of photon production 

If we consider the case of hydrogen, then there is a principle dif- 
ference in the efficiency of (low frequency) photon production if 
the feedback process occurs in the H i Lyman continuum or the Ly- 
man series. The efficiency of photon production will be larger if the 
feedback photon disappeared in the Lyman continuum, while for 
feedback in the Lyman series the total number of (low frequency) 
photons is not changing very much. Therefore it is also clear that 
the total efficiency of secondary photon production will depend on 
whether the Hi Lyman continuum is optically thick or thin. A sim- 
ilar statement also applies to the photon production efficiency of 
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Figure 14. Death probabilities in tlie Hi Lyman series. We included 200 
shells into the sums and assumed that the ambient radiation field is given by 
the CMB blackbody spectrum. 




Figure 15. Relative contribution of transitions towards lower levels to the 
total death probabilities in the H i Lyman series. We included 200 shells into 
the sums and assumed that the ambient radiation field is given by the CMB 
blackbody spectrum. 



He I. Furthermore, the same physical argument explains why the 
total number of photons is not affected very much in the case of 
intra-species feedback. 

To understand this point, let us consider a feedback photon 
that dies in the H i Lyman continuum. It liberates an electron which 
then can be captured to some highly excited state, sub-sequently 
emitting several low frequency photons on it way back to the 
ground state. Here the important point is that it in principle can 
reach states with large angular momentum quantum number. These 
states can only depopulate via many transitions with small hn, i.e. 
preferably emitting low frequency photons. 

If on the other hand, a photon is feeding back on one of the H i 
Lyman series transitions, then even for Ly-« with n> 2 the electron 
will very likely remain bound and only reach a few intermediate 
states preferably with small angular momentum quantum number, 
before again ending in one of the Lyman resonances. The low /- 
states will nearly always depopulate via a few transitions with large 
An, i.e. only emitting some photon at high frequencies. 

To support these statement, in Fig.[l4]we give the death prob- 
abilities for several Hi Lyman series transitions. In Fig.ll5landll6l 
we also show the partial contribution to the death probably com- 
ing only from transitions towards lower levels and the continuum, 
respectively. Even though for levels with n > 2 the death probabil- 
ity is very large in comparison to the H i Lyman a line, less than 
15% - 25% of the transitions end with a free electron. In all cases 
about 40% - 50% of transitions go towards lower levels, and the 
remaining ~ 25% - 35% lead to some higher levels. 

As we will see below (Sects. 18. 9t , this aspect of the problem 
changes the total number of secondary low frequency photons that 
are produced due to feedback of Heii photons. 



8.4 Possible speed-up of He hi He ii recombination by He i 

Like in the case of Heii — » He: recombination, the absorption 
of Heii photons by the small amount of Hei already present at 
z ~ 6000 could speed up He hi — > He ii recombination (5000 < z < 
7000). However, there are several differences that make this pro- 
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z 

Figure 16. Relative contribution of continuum transitions to the total death 
probabilities in the H i Lyman series. We included 200 shells into the sums 
and assumed that the ambient radiation field is given by the CMB blackbody 
spectrum. 



cess much less important than in the case of hydrogen. Although 
the dependence of the photon-ionization cross section for neutral 
helium on frequency is not as strong as in the case of hydroge 
jSamson et al.lll994l) . the amount of helium is about ~ 13 times 
smaller than the amount of hydrogen. Also during the He hi — > He h 
recombination He i is in addition destroyed, because He h becomes 



In the energy range 



24.587 eV < hv < 55 eV we found the fit 



cr^^' « 7.4 X 10"'** cm^ [v/v!^"]-'■" to the data of ISamson et"ai] fT9» 



Is 

which is accurate to ±5%. 
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Figure 17. Absoq^tion optical depth in the Hei I'Sq continuum for a 
He II Lyman a photon that was emitted at redshift z and traveled until 
z' = z - Az. At high redshifts (z > 4500) we used the Saha solution for 
the He I 1 ' Sp populatio n and at low redshifts {z < 4500) the Recfast solu- 
tion )SeageretalJl999h . The curve labeled 'Az/z — > lim' denotes the cases 
when the photons reach below the continuum threshold. 
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Figure 18. Absorption optical depth in the Hi Is continuum for a Heh 
Lyman a photon that was emitted at redshift ; and traveled until z' = z- 
Az- At high redshifts (z > 4500) we used the Saha solution for the Hi Is 
population and at low redshifts (z < 4500) the solution obtained with our 
implementation of Recfast. The curve labeled 'Az/z —> lim' denotes the 
cases when the photons reach below the continuum threshold. 



ionized. This effect is also increased because the ratio of the ioniza- 
tion potentials of neutral helium and singly ionized helium is a bit 
larger than for neutral helium and hydrogen. This leads to an addi- 
tional decrease in the amount of neutral helium close to the time of 
He III He ii recombination. 

To understand the importance of the He: I'Sq continuum 
opacity we computed the the absorption optical depth for the He ii 



Lyman a line for different Az/z between the emission redshift z and 
the absorption redshift z' = z - Az < z (see results in Fig.ll7ll. For 
comparison, the Doppler width of the the He ii Lyman a line is of 

[-■1/2 
. From Fig.[T7]we can see that 

inside the He n Lyman a Doppler core the He i 1 ' So continuum ab- 
sorption should become important only at z ~ 3000 - 3250. At that 
time practically all He ii recombined, so that there is no important 
speed-up of He m — > He ii recombination expected. 

Comparing more carefully with the speed-up seen for Hei, 
from Figs. [6] and [3] it is clear that even t^' ~ 10"'' already leads 
to some notable effect. However, in for the Hei I'Sq continuum 
absorption inside the Heii Lyman a Doppler core this is reached 
only at z ~ 4000 - 4300, again well after Heii recombined. In 
addition, Heiii — » Heii recombination already occurs very close 
to the Saha case, so that the shape and width of the spectral lines 
cannot be altered much in addition. We therefore believe that this 
process is negligible in this context. 

Finally, also the small amount of neutral hydrogen present al- 
ready during He iii — > He ii recombination cannot change anything 
about this conclusion, since inside the He n Lyman a Doppler core 
the H I continuum optical depth it completely negligible (cf. Fig.|18l 
curve for Az/z ~ 10"^). 



8.5 Electron scattering during He iii He ii recombination 

Although in terms of the ionization history during hydrogen and 
neutral helium recombination electron scattering is not very impor- 
tant, for the shape, posit ion and width of the Heii distortions its 
effect is already notable l lRubino-Martm et"^l2008l) . This should 
affect the time-dependence of the feedback caused by high fre- 
quency Heii photons on Hei, since, for example, electron recoil 
will lead to some important drift of the lines towards lower frequen- 
cies, increasing the effective rate of redshifting. Also the Doppler 
broadening of the lines will affect the duration of the feedback. Fur- 
thermore, one also expects that electron scattering could lead to an 
significant increase in the Heiii — > Heii recombination rate. For 
the ionization history and computations of the CMB power spectra 
this again will not be important, but in computations of the CMB 
spectral distortions introduced during recombination this process 
will matter. Also one shoul d include the CMB induce d stimulated 
electron scattering process l lChluba & Sunvaevl2008ah for accurate 
computations of the low frequency spectral distortions. 



8.6 Time of He ii Lyman a feedback on He i 

From Fig.[T7]one can understand when most of the Hed Lyman a 
feedback is expected to occur. Since most of these photons were re- 
leased at z ~ 6000, we can see from Fig.[T7]that they will disappear 
in the Hei I'Sq continuum until Az/z ~ 0.3, or z' ~ 4200. This 
is well after the recombination epoch of He ii and still well before 
the normal recombination epoch of He i. Therefore the feedback is 
most strong in the pre-recombinational epoch of neutral helium. 

We computed the strength of this feedback in more detail us- 
ing the result for the Heii Lyman a line from our 100 shell com- 
putation. The resulting correction to the Hei I'So photoionization 
rate is shown in Fig. [19] As one can see the feedback is largest at 
z ~ 4280. The width of the feedback curve is about 10% - 15%, so 
that one can expect some rather narrow pre-recombinational, feed- 
back induced lines in the He i recombination spectrum. 

Looking again at Fig. [17] one can also see that practically 
all He II Lyman a photons will never reach frequencies below the 
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Figure 19. Feedback induced correction to the Hei I'So photoionization 
rate. We only included tlie effect of the He ii Lyman a line in our computa- 
tions using the pre-computed results from our 100 shell calculations. 




Figure 20. Absorption optical depth in the Hi Is continuum for a 
Hei 2^Pi - I'So photon that was emitted at redshift z and traveled until 
z' = z — Ai. At high redshifts (z > 4500) we used the Saha solution for the 
H I 1 s population and at low redshifts {z < 4500) the solution obtained with 
our implementation of Recfast. The curve labeled 'Az/z — * lim' denotes 
the cases when the photons reach below the continuum threshold. 



He I 1 ' So continuum threshold. This is because practically no pho- 
tons are released at z > 7000, and there the maximal optical depth 
(corresponding to the case Az/z — > lim) exceeds unity already. 

One should also mention that the a small fraction of He ii Ly- 
man a photons is expected to directly disappear in the H i Lyman 
continuum (see Fig.|18K but this should have a negligible effect on 
the recombination spectrum. 



8.7 Reprocessing of pre-recombinational He i photons 

As explained above (Sect. [8^ . every Heii photon that was repro- 
cessed by He i will eventually be replaced by some photon in res- 
onances connecting to the He i 1 ' So state. Since at z > 2500 the 
Hi continuum opacity between these Hei resonances is not im- 
portant (cf. Fig[3}, it is clear that for these secondary Hei pho- 
tons one only has to look at the possible feedback process from 
the Hei 2^Pi - I'So intercombination line. 

In Fig.|20]we present the optical depth in the Hi Lyman con- 
tinuum for photon emitted at the different redshift in the He i 2^Pi - 
I'So line. Again we consider several values for the considered red- 
shift interval that the photon traveled. It is clear that photons emit- 
ted at z < 3700 - 4000 will all be absorbed in the H i Lyman con- 
tinuum before they can reach below the threshold frequency (cf. 
curve labeled 'Az/z — > lim'). However, most of the secondary Hei 
photons emitted at z ^ 4000 will in fact survive the absorption in 
the H I Lyman continuum and instead die in the H i Lyman series 
at lower redshift. This therefore leads to a very extended feedback 
processes of secondary He i photons on hydrogen. As we will see 
below (Sect. [8^ the associated spectral distortion at high frequen- 
cies indeed is rather broad, reflecting this fact. 

8.8 Changes in the recombination spectrum due to He n 
Lyman a feedback 

In this Section we want to demonstrate the effect of Hed feedback 
on the cosmological recombination spectrum. Taking the whole 
list of processes mentioned above into account is far beyond the 
scope of this paper, however, it is rather straightforward to account 
for the feedback effect from the Heii Lyman a line, using pre- 
computed results for the changes in the He i 1 ' So photoionization 
rate (see Fig. [19] and Sect. l8.6l >. The main results of these compu- 
tation are shown in Fig. 12 11241 The most important aspect is that 
there are several additional, rather narrow spectral features in the 
He I bound-bound spectrum (cf. Fig. I23t . which locally increase 
the contribution coming from helium significantly (see Fig. 124b . 
This may render a determination of the primordial helium abun- 
dance using the cosmological recombination spec trum (e.g. see 
IChluba & Sunvaevll2008b| - |Sunvaev & Chluball200^ a bit easier. 

8. 8. 1 Changes in the Hi recombination spectrum due to He ii 
Lyman a feedback 

As in the case of Hei feedback on hydrogen, one can see (see 
Fig. 1211 ) that the feedback of quanta related to He ii leads to the 
production of several additional photons in the recombination spec- 
trum of hydrogen. In particular, the pre-recombinational feature at 
high frequencies (see Fig. [21] middle panel) changes significantly, 
while a low frequencies, the additions are generally smaller (see 
Fig. [2T] upper panel). Only at frequencies around v ~ 200 GHz, 
one can see significant changes. As explained above (Sect. 18. 6I ). 
He n Lyman a photons do not directly feed back on hydrogen, since 
these quanta are all absorbed in the He 1 1 ' So continuum before. 
Furthermore, most of the secondary Hei high frequency photons 
feed back on the Hi Lyman series (Sect. lSJt . so that, in agreement 
with our computations, the low frequency photon production effi- 
ciency is smaller (see Sect. l8.3t . The free-bound spectral distortion 
(see Fig.[2T] lower panel) like in the case of He i feedback (e.g. see 
Fig. 1 12b again shows practically no additional features, but only a 
small additional increase in the total amplitude. 

In Fig. [22] as an example, we also demonstrate the changes 
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Figure 21. Contributions to the cosmological recombination spectrum from 
tiydrogen. We sliow the low frequency Hi bound-bound spectrum (upper 
panel), the high frequency Hi bound-bound spectrum (middle panel), and 
the Hi free-bound spectrum (lower panel). We included 20 shells of hydro- 
gen and 20 shells of helium in our computations and allowed for helium 
feedback from the first ten shells. The curves marked with 'I^CS 2008 ' 
where computed using our reference model iRubiiio-Marti'n et alj|2008l) . 
Those marked as 'normal recombination spectrum' do not include the ef- 
fects of the hydrogen continuum opacity on the dynamics of helium recom- 
bination. For the cases with He ii feedback we only included the effect of 
the He n Lyman a line and allowed for escape in the H i and He i continua. 
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Figure 22. Feedback induced correction to the Hi Lyman a line (upper 
panel) and Hei 2'Pi-l 'So line (lower panel). We only included the effect of 
the He n Lyman a line in our computations. We show the spectral distortion 
Aly as a function of the emission redshift, which can be computed with the 
relation v(; = 0) = v,7[l -l-Zem]- In both cases we subtracted the phase space 
density on the fai' blue side of the resonances. Furthermore, we show the 
cases with and without allowing photon escape in the H i and He i continua. 



in the H i Lyman a line as a function of redshift. In addition to the 
cases shown in the previous figures, we also ran a case for which we 
did not include the escape of photons in the H i Lyman and He 1 1 ' So 
continuuin. Although for the distortions from hydrogen this does 
not make a large difference, for the distortions from helium it is 
very important (see next Section). 



8.8.2 Changes in the He: recombination spectrum due to Heh 
Lyman a feedback 

In Fig.|23]we show the Hei bound-bound recombination spectrum 
for different cases of the feedback process. One can clearly see that 
the feedback of He ii Lyman a photons affects the He i emission 
very strongly, leading to several new features at practically all fre- 
quencies. It is also clear that the escape in the Hei I'So continuum 
is very important for the correct shape and strength of the pre- 
recombinational features. Without the inclusion of Hei 1 'So con- 
tinuum escape the pre-recombinational features are overestimated 
by about a factor of 2. This again is due to the fact that every feed- 
back photon that can escape the Hei 1 'So continuum will subse- 
quently feed back on bound-bound transitions connecting to the 
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Figure 23. Contributions to the cosmological recombination spectrum from 
neutral helium. We show He i bound-bound spectrum. We included 20 shells 
of hydrogen and 20 shells of helium in our computations and allowed for 
helium feedback from the first ten shells. For the cases with Heh feedback 
we only included the effect of the He it Lyman a line. 

ground state of helium, for which the low frequency photon pro- 
duction rate is expected to be smaller (see explanation in Sect. l8.3t . 

As an example, in Fig. [22] we give the distortion due to the 
He 1 2' Pi - I'So resonance only. One can see that the maximum 
of the pre-recombinational feedback is reached at z ~ 4200, cor- 
responding to the maximum in the correction to the He: I'Sq 
photoionization rate (see Fig. |19ll. Figure [22] also shows that the 
pre-recombinational feature from a given resonance is expected to 
have a larger amplitude than the recombinational feature itself, even 
though the total emission remains comparable. This is because the 
feedback process occurs over a rather short time interval, owing to 
the narrowness of the FIe ii Lyman a line. As explained in Sect. |8.5l 
the inclusion of electron scattering into the problem may alter this 
aspect of the problem. In addition, the number of ionizing Heii 
photons per He i atom is comparable to unity, so that for full feed- 
back in the He 1 1 ' So continuum the number of He i quanta would 
practically be doubled. However, as mentioned above the efficiency 
of low frequency photon production is decreased due to the partial 
escape in the He 1 1 ' So continuum. 

8.8.3 Changes in the total H i and He i bound-bound 
recombination spectrum 

Finally in Fig. [24] we present the sum of the Hi and Hei bound- 
bound recombination spectrum. We only included the contributions 
due to dipole transitions, but neglected those from the free-bound 
and two-photon continua, since they are not leading to some vari- 
able signal and hence are not as interesting from an observational 
point of view. One can clearly see that the feedback of helium pho- 
tons leads to non-trival modifications in the recombinational radi- 
ation, with the strongest helium-related features at v ~ 200 GHz, 
~ 80 GHz, ~ 35 GHz, and ~ 10 GHz. However, details in the spec- 
trum appear at practically all frequencies, and in particular the vari- 
ability of the spectral distortion is altered, with significant changes 
in the amplitude, positions of peaks and new features appearing. 
All these traces in the recombination spectrum are only due to the 
presence of helium in the early Universe, and therefore may render 
it possible to determine the primordial helium abundance using the 
spectrum of the CMB. Most interestingly, the amount of helium- 
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Figure 24. Total bound-bound recombination spectrum from H i and He i 
at low frequencies. We only included the spectral distortions due to dipole 
transitions. 

related photons is increased when accounting for all the feedback 
processes. This should make it slightly easier to use the spectrum 
of the CMB for measurements of the helium abundance. 

8.9 Estimates for the total number of feedback induced 
photons from helium and hydrogen 

Using the spectral distortions discussed in the previous section, one 
can now compute the total number of photons emitted by the dif- 
ferent atoms. In Table[3]we summarize the results of these compu- 
tations. In total one finds about 3.2 additional photons per helium 
atom when including all the feedback processes. This corresponds 
to an increase in the total number of helium-related photons by 
~ 37% in comparison to the standard calculation which assumes 
independent recombination histories. I n comparison to the refer- 
ence model jRubifio-Martm et"al]|2008h , the net gain is about 2.38 
additional photons per helium atom, so that the feedback induced 
number of photons increased by a factor of 3.2/0.77 ~ 4.2. 

Furthermore, per helium atom one has a total of ~ 1 1.5 y re- 
lated to the presence of helium in the Universe, or about 0.91 y per 
hydrogen atom, resulting in a total contribution of helium-related 
photons to the recombination spectrum of ~ 17%. In the standard 
calculation which assumes independent recombination histories he- 
lium contributes about 13%. However, here one should mention 
that this average increase in the helium-related number of photons 
can be exceeded several times at particular frequencies (see Fig.l24l 
and Sect. [8X3] l. 

Since for our computations in total 0.9 -I- 0.55 ~ 1.45 quanta 
per helium atom were re-processed, the above results also imply 
that per feedback photon about 3.2/1.45 ~ 2.2 additional photons 
were produced. The production only due to primary He i photons 
IS abou0 1.9/0.9 ~ 2.1 additional photons, while the production 
due to Heii photons is about 1.3/0.55 ~ 2.3 additional photons per 
helium nucleus. Although from naive estimates due to the double 

Note that here we also included the small reduction (by ~ 0.3 y per 
hehum nucleus) of the total number of He i photons, which is due to the 
changes in the dynamics of helium recombination caused by the H i contin- 
uum opacity. As explained in Sect. |7.1.1| per primary He i feedback photon 
Hi produces ~ 2.6 extra photons per helium nucleus. 
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Table 3. Total number (density) of photons, Ny, produced by dilferent 
species. Here N}i and Nne denote the total number (density) of hydrogen 
and helium nuclei. Contributions marked with 'rec' are from the normal re- 
combinational epoch. Those mai'ked wit h 'pre-rec F are the addition al pho- 
tons computed for our reference model iRubiiio-Marti'n et alJIlOOSh . while 
those marked with 'pre-rec IF give the additional photons when accounting 
for the full Hei feedback. Those marked with 'pre-rec 111' denote photons 
that aie produced when including He n Lyman a feedback. We also give the 
numbers for the whole bound-bound and free-bound spectmm, as well as 
the total difference (bound-bound + free-bound) in the number of photons. 
We included 20 shells for H i and 20 shells for He i in the computations. For 
Heh we used the results for 100 shells. 



Atom 






ti I (rec ) 


4 51 




H I (pre-rec 1) 


n 1 9 


1 .J 


H I (pre-rec II) 


0.12 


1.5 


H I (pre-rec III) 


0.1 


1.3* 


H I (all) 


4.85 




He I (rec) 


0.29 


3.69 


He I (pre-rec 1) 


-0.06 


-0.73+ 


He I (pre-rec II) 


-0.04 


-0.45 


He I (pre-rec HI) 


0.04 


0.53* 


He I (all) 


0.24 


3.04 


Heh (rec) 


0.37 


4.69 


He n (pre-rec III) 


-0.04 


-0.55* 


Heh (all) 


0.33 


4.14 


total number (rec) 


5.17 




total difference (pre-rec I) 


0.06 


0.77 


total difference (pre-rec II) 


0.08 


1.1 


total difference (pre-rec III) 


0.10 


1.28 


total number (all) 


5.42 




total difference (pre-rec I-III) 


0.25 


3.2 



* This number includes all surviving high frequency photons from He i and 
those from the H i Lyman continuum. 

* About -0.45 is due to the removal of high frequency photons in the H i 
Lyman continuum. The rest is due to changes in the low frequency photon 
production as a result of the speed-up of He i recombination. 

* This number is only an estimate, since some of the 0.55 ionizing He i pho- 
tons directly escaped in the He 1 1 ' So continuum. We accounted for these 
photons in the computations, but did not check their number exactly. The 
difference in the Hei bound-bound spectmm was ~ 0.25 y per helium nu- 
cleus. To this we added 0.55 X 50% ~ 0.28 of high-n free bound transitions 
resulting in the quoted number 

* We assumed that only the He ii Lyman a line is re-processed. 



re-processing of He ii photon one would expect the efficiency of 
photon production to be larger in the latter case, due to the fact 
that the hydrogen and Hei Is-continua are not completely optically 
thick (see Sect.[8j3]l the numbers remain comparable. 

With these values one can also provide an optimistic upper 
limit on the total amount of helium-related photons. Assuming that 
every high frequency photons emitted by Heii lead to ~ 2.3 ad- 
ditional photons, using the numbers given in Sect. 18.11 we expect 
2.3 X (0.98 -I- 0.5) ~ 3.4 additional photons due to Heu feed- 
back alone. This is another factor of ~ 2.6 times more than we 
already obtained here. Also adding those photons that may be 
(see Sect. |7.lT| and [7!2l > produced due to primary Hei feedback 
(~ 2.6 - 0.3 ~ 2.3 y per helium nucleus) one can therefore expect 
a total addition of ~ 5.7 y per helium atom, when including all the 



feedbacks. This implies a total of ~ 14 y related to the presence of 
helium in the Universe, or about 1.1 y per hydrogen atom, resulting 
in a total contribution of helium-related photons to the recombi- 
nation spectrum of ~ 20%. In comparison to the computation that 
does not include any of the feedback processes one therefore could 
expect a factor of 14/8.4 ~ 1.7 times more helium-related photons. 



9 SUMMARY AND CONCLUSIONS 

In this paper we considered the re-processing of high frequency 
photons emitted by helium during the epoch of cosmological re- 
combination. We investigated both the possible changes in the cos- 
mological recombination spectrum and the cosmological ionization 
history of the Universe, taking several previously neglected physi- 
cal processes into account. In particular we studied 

(i) the feedback of He i photons on He i recombination 

(ii) the feedback of primary He i photons on H i recombination 

(iii) the feedback of primary He ii photons on He i recombination 

(iv) the feedback of secondary Hei photons on Hi recombination. 

The proce s s (i) was already considered earlier 
ISwitzer & Hirat j ( l2008j) . but the processes (ii)-(iv) have not 
been studied at full depth so far. We find that process (i) leads to 
no important change in the cosmological recombination spectrum, 
but a small modification in the cosmological ionization history by 
AA'e/A'c ~ +0.17% at ; ~ 2300. At low redshift our correction 
seems to be smaller (s ee Fig. Hot than the one presented in 
ISwitzer & Hirat j ( l2008hl) . However, it is clear that the difference 
will not be very important for the analysis of future CMB data. 
Nevertheless, one should check the reason for this difference in 
a detailed code comparison, which is currently under discussion 
among different groups. 

Furthermore, we find that the processes (ii)-(iv) affect the ion- 
ization history at a negligible level, but could increase the total he- 
lium induced contribution to the cosmological recombination spec- 
trum by ~ 40% - 70% in comparison to the standard calculation 
which assumes independent recombination histories with no feed- 
back between the different atomic species (see Sect. |8.9l >. 

The first result is related to the fact that the feedback processes 
(ii)-(iv) all occur in the pre-recombinational epochs of the consid- 
ered atomic species, where they are still very close to full equi- 
librium with the free electrons and protons, so that perturbations 
in the ionization state are restored very fast. On the other hand, 
the increase in the total number of helium-related photons occurs 
because the recombination of electrons that were liberated by the 
feedback process can be captured into highly excited atomic lev- 
els with the possibility to release several photons per one ionizing 
photon in the subsequent cascade to the ground state. Here in par- 
ticular the double re-processing of He ii photons (see Sect. 18. 2t is 
interesting, for which primary He i photons in principle feed back 
on both He i and at later times on H i. Therefore one naively expects 
the feedback induced number of photons from He ii to be twice the 
one coming from Hei feedback. However, for Heii feedback the 
efficiency of low frequency photon production is not as large (see 
Sect. 18. 3t , because the hydrogen and neutral helium is continua 
are not completely optically thick at the time of feedback. There- 
fore the final number of produced photons per ionizing He ii photon 
remains comparable to the one for Hei feedback (see Sect. 18. 9b . 

As a result of all the feedback processes included here a to- 
tal of ~ 12 y per helium nucleus are produced during cosmological 
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recombination, implying a total contribution of ~ 17% to the cos- 
mological recombination spectrum. In total about ~ 3.2 additional 
helium-related photons per helium nucleus were released in com- 
parison the standard calculation which neglects feedback. This is 
an important addition to cosmological recombination spectrum and 
may render it easier to determine the primordial abundance of he- 
lium using differenti al measurements of the C MB energy spectrum 
(e.g. as discussed in ISunvaev & Chlubj|2009l) . Here most interest- 
ingly, the feedback of Heh photons on Hei leads to the appear- 
ance of several additional rather narrow spectral features in the He i 
recombination spectrum at low frequencies (cf. Fig|24t. The con- 
tribution of these features to the spectral distortion at some given 
frequency can exceed the average level of 17% several times. In 
particular the bands around v ~ 10 GHz, ~ 35 GHz, ~ 80 GHz, 
and ~ 200 GHz seem to be affected strongly (see Fig|24b. In ad- 
dition, one should point out that today these photons are the only 
witnesses of the feedback process at high redshift, and observing 
them in principle offers a way to check our understanding of the 
recombination physics. 

Since for the feedback of He ii photons here we only accounted 
for those from the Heii Lyman a resonance, which is only ~ 1/3 
of the total number of primary Heii photons potentially available 
for feedback (see Sect. 18.11 for more detailed estimates), one can 
still expect the total number of helium-related photons to increase 
by another ~ 10% - 20% to a total of ~ 14 y per helium nucleus 
(see Sect. |8.9| for estimates). However, we leave a detailed investi- 
gation of the He ii feedback problem including the two-photon con- 
tinuum, Balmer continuum, and the effect of electron scattering for 
some future work (see Sect. [8] for a more detailed overview of the 
possible processes related to this problem). Also in the Hei feed- 
back problem so far we neglected the effect of the Hei2'So - I'So 
two-photon continuum and the ionizing photons from levels with 
n > 10. However, their effect on the cosmological recombination 
spectrum is expected to be smaller (see Sect. |7.2| l. 

Finally, we also refined our previous compu tations of He ii 
He: recombination jRubifio-Martm et"al]|2008l ) checking the ef- 
fects of (A) Hei ;j'D2 - I'So electric quadrupole transitions; (B) 
higher (n > 2) He: n^V\ - I'So intercombination transitions; 
(C) the speed-up of the higher (w > 2) He: ii'Pi - I'Sq se- 
ries due to the absorption of helium photons by hydrogen; (D) 
possible time-dependent corrections to the Hei 2'Pi - I'So and 
He: 2'Pi - I'So channels; (E) the effect of the thermodynamic fac- 
tor on the He: 2' Pi - I'So and He: 2^Pi - I'So transitions. These 
processes affect the ionization history but do not introduce any sig- 
nificant changes in the recombination spectrum. 

For the processes (A)-(C) we find very good agreement with 
the results of lSwitzer & Hiratj l l2008d) . and conclude that they can 
be neglected in the analysis of future CMB data. Although pro- 
cesses similar to (D) and (E) le ad to significant correction s during 
the epoch of H i recombination jChluba & Sunvaevll2009dl lbl). here 
we show that for the recombination of helium they are of minor 
importance, mainly because the effect of the Hi continuum opac- 
ity on the Hei 2'Pi - I'Sq and Hei 2'Pi - I'So channels is much 
more crucial. Therefore, they again can be neglected for theoretical 
computations of the CMB power spectra. 
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APPENDIX A: ATOMIC MODEL FOR HELIUM 

We use the sam e ato mic model for helium as in 
iRubifio-Martm et al] | |2008|) . which is mostl y based on en- 
ergy levels and transition rates given iDrake & MortonI 
( |2007|) . However, here we also took additional He: «^Pi - I'So 
intercombination transitions with 3 < ;i < 10 into accou nt. For the 
He I 3 ''Pi - 1 ' So transition we used the value A = 56.1 s"' jLaughliiil 
Il97^ and like in Swi tzex& Hirata ( 2008a) we extrapolate to larger 
n assuming A oc I/n^. Furthermore, we also added electric 
quadrupole transitions (E2) for the sequence Hei jj'Dt - I'So 
with n < 10 using the absorption oscillator strengtlP^ given by 
ICann & Thakkad ( l2002h f or n < 6 an d extrapolating these values 
assuming / oc 1 /n^ (like in ISwitzer & Hiiata.,2008^) . 



APPENDIX B: ESCAPE IN THE CONTINUUM 

Integrating Eq. l|6j over dv df2 one can directly write 



r 1 dNy 

J c "dT 



dvdn = N^Rll - + ARi^l (BI) 



Here one has 



■■ 47t[N[ 
: 4n[N[ 



:„,-<'] j o-i,,(v)G;;dy 



(B2) 



Since in the Wien tail of the CMB one has f ^^^^ dv = 
crwr(v) Hir dv = — one can introduce the normalized profile 

V^y^^ z) = /%" ™ ^ write 



«■ Ai?„, = %[yv:„, - <] r ^,(v, z) n Gl dv. 

= [N, <^ - A'|1'<J J Vc(y, z) /; Gl dv, (B3) 



where we used the fact that ^ fi^{T)lf^^^lN^}^ = if^^^. Therefore 
one finally has 



J c df 

with /-^^ = X" Vc(v, z)[l - /; Gl] dv. 



dv dQ. = \N, N, Rl\^ - N'^: R'^ X Pli, (B4) 
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i200^ and references therein. 

'' These are related to the quadrupole transition rate by (e.g. ISobelmanl 
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E2 _ Sn^e-' 
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Bl Quasi-stationary solution witliout asymmetry between 
continuum emission and absorption profile 

If we assume quasi-stationary conditions for the continuum and ne- 
glect tlie factor l/f^ in the definition of &l, Eq. l |4bt . then one di- 
rectly has G" ''' = 1 - e-'-'^*""-*. Using Eq. l|7bj and Eq. ^ with 
the substitution i' = v -j-j^ one can furthermore write 



V 



(B5a) 
(B5b) 



Here t^'"^ = -tt^ « 5 — 1^ r-^ and = ii>Av )dv . 

Also in Eq. dBSbt we have neglected the factors y and /^(z). For 
the escape probably from the main resonances in helium and hy- 
drogen, this would be the st andai'd procedure to obtain the Sobolev 
escape probability (e.g. see IChluba & Sunvaevll2009d^ . Then with 
^Pc(y^ /v dv after also setting ~ \ one finally 



DC 

has 



rjC.q.s.noasym 
-^cff 



(B6) 



This expression has the same form as the standard Sobolev escape 
probability. However, here we have neglected the fact that the emis- 
sion and absorption profiles in the hydrogen Lyman continuum dif- 
fer by fl, which makes the absorption profile iff'^iV) = fyifdv) ~ 
CiscCv) very broad. This is the reason why the expression for the es- 
cape probability in the continuum scales like P^T ~ 1/ [1 + Tf^] in- 
stead, as already obtained in Ichluba & Sunyaev] l l2007h . Below we 
give a short derivation of this expression in the formulation given 
here. 



B2 Quasi-stationary solution including the asymmetry 
between continuum emission and absorption profile 

Here we give the derivation for the escape probability in the Hi 
Lyman continuum assuming quasi-stationarity but taking into ac- 
count that the continuum emission and absorption profiles differ by 
t I. From the definition Eq. jlhj and with v = v it directly 
follows 



T',t,(v,z',z) ~ To 



r 



O-lsc(v) V 



■ To ■ 



(B7) 



with To 
has 



II T — 



Also from Eq. ( I4bb with 0J^ = 1 one directly 

(z) cri,,(v') 



To + 



exp 

hv^ 



-To ■ 



dv' 
v' 
dv' 
v' 



1 Tf'-Xv) 

^ TTtP(v)' 



(B8) 



with Tf^(v) = ""'"']J'"'' Inserting this into the expression for 
the escape probability P^^^ = (/)c(v, z)[l - f^C^^Av, and ne- 
glecting the slow (power-law) scaling of t^^'^(v) » t^'"^(vc) with fre- 
quency one then obtains » 1/[1 -1- rf^]. 

The main difference to the approximation Eq. l lB6t is in the 
scaling for small rf^. Here one has a 1 
Eq. ^ one finds p^;^^^ ™-^- « 1 - 
t^*"^ » 1 both approximations agree. 



■^■^ while from 
However, in the limit 
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